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Abstract 

We discuss the use of field theory for the exact determination of universal properties in two- 
dimensional statistical mechanics. After a compact derivation of critical exponents of main 
universality classes, we turn to the off-critical case, considering systems both on the whole 
plane and in presence of boundaries. The topics we discuss include magnetism, percolation, 
phase separation, interfaces, wetting. 
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1 Introduction 

Statistical systems close to a second order phase transition point exhibit properties which do 
not depend on details of the microscopic interaction but only on global features such as internal 
symmetries and space dimensionality. These properties are called universal, and systems sharing 
them are said to belong to the same universality class. Renormalization theory has progressively 
clarified how the divergence of the correlation length as the critical point is approached allows 
for the emergence of universality, and how field theory actually is the theory of universality 
classes. 

While these ideas apply in any dimension d > 2, to the point that d serves as an expansion 
parameter in the Wilson-Fisher approach to renormalization, different dimensionalities exhibit 
specific features which are intrinsically non-perturbative and determine distinctive qualitative 
properties of universal behavior. The two-dimensional case, in particular, is characterized by 
features such as absence of spontaneous breaking of continuous symmetries, fermionization (or, 
conversely, bosonization), correspondence between interfaces and particle trajectories, infinite 
dimensional conformal symmetry at criticality. The last property is in turn at the origin of a 
further, remarkable peculiarity of the two-dimensional case, i.e. the existence of exact solutions, 
both at criticality and near criticality, and for all universality classes. 

The aim of this article is to illustrate how field theory leads to the exact description of 
universality classes of classical equilibrium statistical mechanics in two dimensions. The pre¬ 
sentation focuses on basic ideas and examples, with a minimal amount of technical aspects. It 
relies on several recent results, and is original for the spectrum of the topics discussed, ranging 
from percolation to wetting, and the way they are cast into a unified theoretical framework. 

The discussion is structured according to the following path. In the next section we recall 
generalities about universality, field theory and its particle description. In section [3] we give a 
concise overview of two-dimensional conformal field theory suitable for our subsequent discussion; 
in particular, we do not focus on minimal models, leaving room for cases, like percolation, in 
which the bulk correlation functions of the order parameter do not satisfy known differential 
equations. In section [5] we show how the critical indices of main universality classes, including 
percolation and self-avoiding walks, can be identified directly in field theory exploiting the 
insight coming from the particle description. In section [5] we recall how, under conditions 
which turn out to account for the main interesting cases, exact solvability extends to the off- 
critical regime, and how solutions originally obtained in the scattering framework lead to the 
computation of universal quantities such as combinations of critical amplitudes. The extension 
of the formalism to crossover phenomena is illustrated in section [6j In section [7] we turn to 
systems with boundaries, showing how general low-energy properties of two-dimensional field 
theory yield the exact characterization of phase separation and of the interfacial region, an 
analysis extended in section [8] to describe the interaction of an interface with the boundary and 


2 



the effect of boundary geometry. In section [9] we show how the particle formalism provides exact 
asymptotic results for percolation on the rectangle away from criticality, where the methods of 
boundary conformal field theory do not apply. The last section is devoted to few final remarks. 


2 General notions 


2.1 Universality 


In the framework of classical equilibrium statistical mechanics [I] a system is specified by the 
Hamiltonian whose value is determined by the configuration of the system. The expectation 
value of an observable O is the statistical average over configurations 


(o)=| 5: 

configurations 

where T > 0 is the temperatur^ and the normalization factor (partition function) 


( 1 ) 


Z = ^ e-«/^ (2) 

configurations 

ensures that the average of the identity is 1. We are interested in systems which in the limit of 
infinitely many degrees of freedom undergo a phase transition for some critical value U of the 
temperature. In more than one infinite dimension, this is signaled by an order observable whose 
expectation value (order parameter) vanishes above U and is a function of T below T^. The 
transition is said to be of the hrst order if the order parameter has a discontinuity at U, and of 
the second order (or, more generally, continuous) otherwise. 

Phase transitions are normally associated to spontaneous symmetry breaking. A group G, 
mapping configurations into configurations, is a symmetry of the system if T-L is left invariant by 
the action of G. Configurations of minimal energy (ground states) are mapped into each other 
by G, and below U a G-invariant system chooses the phase dominated by a specific ground state 
as T —)■ 0. 

As r —>■ Tc, if the transition is continuous, universal critical properties emerge, i.e. quan¬ 
titative properties common to systems having the same dimensionality and symmetry group. 
Such systems are said to belong to the same universality class. The Ising model is the basic 
representative of the simplest universality class, that associated to the spontaneous breaking of 
the group G = Z 2 . The model is characterized by the Hamiltonian 


Rising — J ^ ^ j ) — dll ) (3) 

(*d> 

where a “spin” variable ai is assigned to the i-th site of a regular lattice and the sum is restricted 
to the pairs of nearest neighboring sites. Unless otherwise stated, we refer to the ferromagnetic 
case J > 0, for which the model possesses two ground states (spins all plus or all minus); the 

^We adopt units in which ks ~ 1- 
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transformation corresponds to the reversal of all spins, and leaves ([3|) invariant. The order 
parameter (cjj) is proportional to (T^ — T)^ as T —)■ T“, with a critical exponent j3 which is 
universal. This can be contrasted with the value of Tc, which changes, for example, when the 
structure of the lattice changes, or next to nearest neighbors interactions are included. The 
spin-spin correlation function (crjcij) decays as as the distance \i — j\ between the two 

spins goes to infinity; such a decay provides the definition of the correlation length ^ that we 
use throughout this article. Continuous phase transitions are characterized by the fact that the 
correlation length diverges as ^ ~ ~ Tc\-’^ when T ^ T^, so that exponential decay of 

correlation functions is replaced at Tc by power law decay. It is the divergence of the correlation 
length which makes microscopic details irrelevant close to criticality and leads to the emergence 
of universal behavior. The susceptibility y = {N —)• oo is the number of sites) 

diverges as 

x~r±|T-rei-^ r^T±. (4) 

While the exponents n and 7 are universal, the critical amplitudes and r± are not; 
and r+/r_, however, are universal, since non-universal ’’metric factors” cancel in the ratio. 
The (;-state Potts model, defined by the lattice Hamiltonian mm 


^Potts - J ^ ^ 

(hi> 


Si = 1,2, 


( 5 ) 


generalizes the Ising model to the case in which the site variable takes q values, to which we will 
often refer as q “colors”. The Hamiltonian is invariant under the action of the group G = Sq of 
permutations of the colors, and the model possesses q ferromagnetic ground states in which all 
sites have the same color. The expectation value of the spin variable 


(Ti 


— .ry 


a = l,...,q.. 


( 6 ) 


vanishes as long as the symmetry is unbroken and provides the order parameter for the transition. 
A particularly interesting feature of the g-state Potts model is that the partition function admits 
the graph expansion [4] 


^Potts = (X - pf^q^^ , (7) 

{si} G 

where G is a graph obtained drawing bonds on the edges of the lattice, Nb is the number of 
bonds in G, Nb the number of edges without a bond, and p = 1 — a set of connected 

bonds is called a cluster, and Nc is the number of clusters in G, with the convention that a 
site not touched by any bond is also counted as a cluster. The graph expansion makes sense 
of the Potts partition function for continuous values of q. In particular, for g —>• 1 the weight 
— p'^^b Qf configuration coincides with that of the bond percolation problem urn , in 
which edges are occupied with probability p and empty with probability 1 — p. Percolation is 
characterized by the existence of a critical threshold pc such that for p > pc the probability of 
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finding an infinite cluster is larger than zero. Equation ([7|) relates such a “geometric” phase 
transition to the thermodynamic phase transition of the ^ ^ 1 Potts model. 

The basic model exhibiting a continuous symmetry is the vector model defined by the Hamil¬ 
tonian 

^vector — J ^ ^ ' Sj ) (8) 

(hi) 

where the spin variable Sj is a n-component unit vector. The symmetry group is G = 0(n) and 
the case n = 2 is usually called XY model, while the Ising model is recovered for n = 1. The 
partition function of the vector model also admits a graph expansion in which n enters as a 
parameter that can be taken continuous. More precisely, the mapping involves a Hamiltonian 
slightly different from Q (see e.g. 0), but belonging to the same universality class. The graphs 
are configurations of closed paths (loops), in which each loop contributes a factor of n to the 
weight. The limit n —)■ 0 is dominated by the configurations with a single loop and describes 
the self-avoiding polymer problem [S]; we will see in section [THl how self-avoidance emerges as a 
property of the universality class. 


2.2 Fields 


The independence of universal properties on microscopic details suggests that universality is 
most naturally captured within a continuous description. Field theory provides such a continuous 
framework for the study of universality classes (see e.g. [7] and references therein). For a d- 
dimensional system we denote by x = (xi,..., Xd) a point in the Euclidean space and replace 
site-dependent lattice observables by x-dependent fields. In particular, the lattice spin variable 
is replaced by a spin (or order) field, <t(x) € R for the Ising model. Similarly, to the energy 
density UiUj entering ([3|) (the prime indicates summation over nearest neighbors of the site 
i) we associate an energy density field e(x). The Euclidean action (or reduced Hamiltonian) 
M is a functional of the fields which specifies the theory and determines the weight of a field 
configuration as where we keep the notation Z for the partition function X^fieid configs 

in the continuum. 

The action A is invariant under the action of a symmetry group G of the system. A further 
characterization can be achieved starting from critical point actions [9]. The divergence of the 
correlation length and the power law decay of correlation functions at a second order phase 
transition point amount to scale invariance of the critical point field theory. Such a theory 
contains no dimensional couplings and is left invariant by scale transformations x ^ ax; it is 
usually called a “fixed point” theory and in this section we will denote fixed point quantities by 
a subscript FP. Scale invariance leads to the two-point function 


($(xi)$(x2))fp 


constant 
\xi — X2|“^^'>’ 


(9) 


for a scaling field <h(x); this equation serves as a definition for the scaling dimension of the 
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field at the given fixed point. The off-critical theory can now be associated to the action 

A = Afp + T J df‘x£{x), (10) 

where r measures the deviation from Tc- Since the action is dimensionless, r has the dimension 
of an inverse length (a mass in our units) to the power d — and provides the dimensional 
coupling which breaks scale invariance. On dimensional grounds we have ^ oc and 

then the critical exponent v = l/{d — X^). 

A G-invariant field theory contains infinitely many helds with growing scaling dimension 
which transform in the same way under the action of the group. In writing (jlOp we omitted 
infinitely many G-invariant fields with scaling dimensions larger than d. Their conjugated cou¬ 
plings have the dimension of a length to positive powers and become negligible when the system 
is observed over distances much larger than lattice spacing, i.e. in the limit in which universal 
properties emerge. In this sense such helds are called “irrelevant”. The action in which irrelevant 
helds are omitted is called the scaling action and characterizes the universality class. The scaling 
action may contain more than one G-invariant held with scaling dimension smaller than d (such 
helds are called “relevant”). This corresponds to theories in which more than one parameter 
needs to be tuned to achieve scale invariance; these theories describe “multicritical” behavior. 
A theory may contain helds with scaling dimension equal to d (“marginal” helds). Marginality 
may be spoiled by logarithmic corrections produced by off-critical interactions. If this is not 
the case, the addition to the hxed point action of such a “truly marginal” held does not break 
scale invariance and generates a line of hxed points. The recognition that the properties of the 
system depend on the scale at which it is observed and the characterization in terms of hxed 
points and relevant and irrelevant helds form the basis of the renormalization group idea 13 E]. 

We saw how the correlation length critical exponent is determined by X^. On the other 
hand, scaling dimensional analysis gives {cr{x)) oc and then /3 = uXa-. Similarly, the 

susceptibility X = f d‘^x(a(x)a(0)} scales as ^ so that 7 = v{d — 2Xfj). Hence the theory 

accounts for the universality of the critical exponents and for the relations among them, and 
reduces their determination to that of the scaling dimensions of the energy density and order 
helds. Also combinations of critical amplitudes like and r+/r_, being dimensionless, are 

universal. It also follows that, while critical exponents are determined by the hxed point theory, 
critical amplitudes need to be computed within the off-critical theory (1101) . 

2.3 Particles 

We are considering the ordinary case of systems that are isotropic in the continuum limit, so 
that the scaling action (1101) is invariant under rotations. Such an Euclidean held theory dehnes, 
when one of the coordinates is made purely imaginary [xd = it), a quantum held theory with 
spatial coordinates xi,... ,Xd-i and time coordinate t, in which relativistic invariance replaces 
d-dimensional rotational invariance. The quantum theory has the same held content as the 
Euclidean one and correlation functions in the two cases are related by analytic continuation from 
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Figure 1: Elastic two-particle process in 1-|-1 dimensions. Time runs upwards. 


real to imaginary time. The quantum theory also admits a description in terms of relativistic 
particles corresponding to the excitations above a minimum energy (vacuum) state; spontaneous 
symmetry breaking corresponds to the presence of degenerate vacua mapped into one another 
by the symmetry m- 

The particle description of a quantum field theory is encoded by the S'-matrix whose 
matrix elements are the probability amplitudes that a particle state at time t = —oo evolves 
into another state at time t = +oo. In general, relativistic scattering conserves total energy and 
momentum, but not the number of particles. We call “elastic” a scattering process in which 
the number of particles and their masses are preserved. Specializing to the case d = 2, to 
which we turn from now on, we depict in Fig. [T] a two-particle elastic process. We consider the 
case in which the particles have the same mass m. The energy e and the momentum p of a 
particle are related by the relativistic dispersion relation e = \/p‘^ + m?. For the process we 
are considering, conservation of energy and momentum imply that the momenta pi and p 2 are 
individually conserved by the scattering. If the particles carry charges, the total charge may 
be redistributed in the scattering, a possibility that we take into account through the particle 
species labels a, 6 , c, d. 

The scattering amplitude of Fig. [1] is relativistically invariant and depends on the single 
relativistic invariant quantity that can be built out of the two energy-momenta; we take this 
invariant in the form of the square of the center of mass energy 

s = (ei 62 )^ - (pi -k P 2 f ■ ( 11 ) 

The amplitude, that we denote as satishes a number of properties m that we now 

state. First of all the amplitude is an element of the unitary S'-matrix. When continued to 
complex values of the variable s, it is an analytic function whose singularities have a physical 
implication. In particular, as a consequence of unitarity, the minimal energy values (thresholds) 
needed for the production of a new final state correspond to branch points of the amplitude; 
these are located at s = (/cm)^ for the production of fc > 1 particles of mass m. A simple pole 
at s = rh^ corresponds instead to a particle of mass m appearing as a stable bound state of the 
two particles in the initial (or final) state of Fig. [TJ 

The unitarity cuts originating from the branch points we just discussed are taken along the 


7 




Figure 2: Analytic structure of the scattering amplitude in the complex s-plane. Unitarity 
(right) and crossing (left) cuts are shown together with two poles. 


positive real axis in the complex s-plane as in Fig. [2l and the physical value of the amplitude 
corresponds to the limit towards the real axis from above + ie) with s G i?, e — 0"'“). 

Such a path lies on the first (called “physical”) sheet of the Riemann surface associated to the 
cut s-plane; additional sheets are accessed sliding into the cuts. 

A further fundamental property of relativistic scattering is crossing symmetry. It states that 
the amplitude for the direct channel process (Fig. [1] with time running upwards) is related by 
analytic continuation to the amplitude for the crossed channel process (Fig. [1] with time running 
from left to right). In the passage to the crossed channel, the particles h and d, which have the 
arrow pointing in the ‘wrong’ direction, are switched into antiparticles b and d, and have their 
energy and momentum reversed (e 2 ,P 2 —& 2 i—P 2 i leading to s —)• 4m^ — s). The crossing 

relation then reads 

Sliis + ie) = 5j^(W -s-ie), ( 12 ) 

for s real. A consequence of (I12p is that an amplitude inherits from the crossed channel branch 
cuts running along the negative real axis, as well as crossing images of bound state poles (Fig. [2]). 

We finally mention the property of real analyticity, which states that the values of the 
amplitude on the upper and lower edge of a cut are related by complex conjugation; this implies 
in particular that the amplitude is real along the uncut portion of the real axis. It also allows 
one to write the unitarity condition as 

^ Sll{s + ie)S^f{s - ie) = Sachd , 4m^ < s < si, (13) 

e,f 

for physical values of s below the first inelastic threshold si. 

Within the particle framework, fields correspond to operators (we denote them by the same 
symbol) which acting on a particle state produce a superposition of particle states. The Hamil¬ 
tonian and momentum operators H and P of the (1 -|- l)-dimensional quantum system act as 
generators of time and space translations and, in Euclidean coordinates xi and X 2 = it, we have 

$(x) = . (14) 

The field/operator is characterized by the form factors 

J^niPl^-'-^Pn) = (0|4>(0)|pi,...p„), (15) 


where we are considering a single particle species in order to simplify the notation, |0) denotes 
the vacuum state (i.e. the state without particles) and |pi,.. .pn) are asymptotic particle states, 











i.e. corresponding to t = —oo (incoming states) or t = +oo (outgoing states) in the scattering 
picture; in these states the particles are widely separated in space and can be treated as non¬ 
interacting. The asymptotic states are eigenstates of H and P with eigenvalues Cj and 

Yl’i=iPii respectively. The two-point function of a scalar field decomposes over the basis of 
asymptotic states as 


($(x)$(0)) = 


(0|^>(x)^>(0)|0) 

dpi 


n=0 


dpn 


27rei 27rer 




(16) 


where we take advantage of rotational invariance to set xi = 0 and adopt the state normalization 
{pi\P 2 ) = 2vrei5(pi — ^ 2 )- For |x| —)■ 00 the correlator decays exponentially, showing that the 
mass m is inversely proportional to the correlation length. If pf" (p) 7 ^ 0 we say that $( 3 :) creates 
the particle. 


3 Critical points 

3.1 Conformal symmetry 

In writing the two-point function ([9]) we assumed that the field <I> is scalar, i.e. invariant under 
rotations. More generally, in two dimensions we will consider scaling fields $( 3 :) with scaling 
dimension which transform as <I’(0) e“*^*“<h(0) under a rotation by an angle a centered 

in the origin; is the Euclidean spin of the field. It is an important property of field theory 
that products of local fields can be expanded onto an infinite-dimensional basis of local fields. 
In a scale invariant two-dimensional theory the need to preserve dimensional and rotational 
properties gives to such an operator product expansion (OPE) the form 

i>j(x)i>j(o) = 

k 

= E $,(0), (17) 

k 

where we introduced the complex coordinates z = xi + 1 x 2 and z = xi — 1 x 2 , which rotate as 
z e*" 2 ;, z as well as the left and right dimensions A<j,, A<j, such that = A$ -|- A$ 

and 5$ = A,j, — A$ (A$. = Xi, ...); the coefficients are called structure constants. Two 
fields <I>j and <I>j are said to be mutually local if the correlation functions (• • • <I>j(x)<I>j(0) • • •) are 
single valued under the continuation 2 ; ^ This is the case if 

Sk Si Sj Ei Z (f^) 

for all contributing to the r.h.s. of (|17l) . 

A field theory contains a symmetric and conserved energy-momentum tensoj^ T^y{x) with 
scaling dimension Xt^^ = d. In two dimensions the three independent components in the basis 

^In real time, energy and momentum are the integrals over xi of Ttt and Ttx^^ respectively. 
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of complex coordinates are T^z = T, Tzz = T and Tzz = Tzz = —with dimensions (A, A) 
equal to (2,0), (0,2) and (1,1), respectively; the conservation equations read 

dT = de, df = de, (19) 

where d = dz and d = dz- In a scale invariant theory. At = Af = 0 means T = T{z), T = T{z), 
i.e. BT = dT = 0, so that the trace = —40 is a constant; X© ^ 0 then implies that the 
energy-momentum tensor is traceless in presence of scale invariance, a property that actually 
holds also in higher dimensions (see e.g. m- 

Within the theories we are interested in for physical applications all fields are mutually local 
with the energy-momentum tensor. We then have the OPE 

OO 

r(z)d>(0)= ^ z-^-^Ln^o), (20) 

n=—cx) 

where the fields Ln^ have dimensions (A, A) = (A$ — n,A<j>). The L^s then act as shift 
operators in the space of fields graded by A (similarly for the Ln one obtains from T). There 
will be fields (p with lowest dimension [primary fields) characterized by the property Ln4> = 0 
for any positive n; the fields Lnfi with n < 0 are examples of the so-called descendants of the 
primary fi', T = L- 2 I is a descendant of the identity. A primary cp and its descendants form 
an operator family that we denote by [p], and the collection of the operator families forms the 
space of fields. 

An additional property of scale invariant field theories is that they are actually invariant 
under the more general group of conformal transformations, i.e. the transformations which 
allow for a point dependent multiplicative change of the infinitesimal distance element (see 
e.g. [7]); in this sense they are called conformal field theories. In two dimensions infinitesimal 
conformal transformations correspond to changes 5z = f[z), 5z = f[z), where / (resp. /) is 
any analytic function of z (resp. z), in such a way that the group of conformal transformations 
in two dimensions has the essential peculiarity of being infinite-dimensional. This turns out to 
imply HSUS] that the operators entering (I20p satisfy the Virasoro algebra 

[Lni Em] — (^ ^)En+m T Y 2 '^)^n,—m 1 (^I) 

where c is a fundamental parameter of the theory known as central charge. One has = A$<I>, 
L_i$ = (9$ and 

T[z)T[0) = ^ + ^T(O) + -^dT{0) + .... (22) 

The central charge of theories which are the sum of sub-systems non-interacting with each other 
is the sum of the central charges of the sub-systems and, more generally, it can be shown that 
c grows with the number of degrees of freedom of the critical point M- Hence, for example, 
Ising, three-state Potts and XY universality classes have to correspond to increasing values of 
the central charge. 

The algebra ()2ip provides a powerful tool for the study of two-dimensional conformal field 
theories, since the operator families correspond to lowest weight representations of the algebra. 
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The dimension (which is also called “conformal” dimension) of the primary cf) is the low¬ 
est weight, and the subspace of descendants with dimension -|- I is spanned by the fields 
... L-jjCf) with 0 < < ... < jj and ~ hence, generically, the dimension of 

such a subspace is given by the number p{l) of partitions of I (the level) into positive integers. 
A very important role in the theory is played by the reducible representations of the algebra 
m, i-e. the representations [(p] which contain another representation whose primary (po is 
a descendant of (p at some level Iq. The irreducible representation obtained factoring out [i;f)o] 
is said to be a representation degenerate at level Iq, and (p a degenerate primary. Clearly, in 
such a degenerate representation the number of independent fields at level I > Iq is smaller than 
p{l). Factoring out [(/iq] amounts to setting to zero (pQ and its descendants, a condition which 
can be shown to lead to linear partial differential equations for multi-point correlation functions 
containing a degenerate primary |12j . 

The basic universality classes of two-dimensional critical behavior correspond to conformal 
theories with central charge c < 1. As pointed out in m, these can be described starting from 
the theory of a free scalar field with action 

A = ^ / . (23) 

The dimensionless character of p{x) leads to the logarithmic correlator 

{(p{x)(p{0)) = — In \x\ = —^(Inz -|- Inz), (24) 

which is consistent with the equation of motion ddp = 0 and the decomposition 

(p(x) = xiz)xiz). (25) 

The components of the energy-momentum tensor can be written in the form 

T{z) =+ iQ , f{z) = -{dp)‘^ + iQd‘^p; (26) 

the presence of the terms containing the parameter Q is allowed by the fact that, due to the 
equation of motion, they do not affect conservation and are total derivatives with respect to xi, 
so that they do not contribute to energy and momentum. Using (|24l) and Wick theorem, we can 
compute T{z)T{w) and, comparing with (1221) . obtain the central charge 

c = 1 - 6Q^ (27) 

Instead of X) which has (x('2)x(0)) oc Inz, proper scaling primary fields of the theory are the 
exponentials 

Vp{z) = , (28) 

whose dimension Ap is obtained considering T{z)Vp{w) and isolating the term ApVp{w){z — w)~‘^ 
(the term with n = 0 in (1201) ): this yields 

A = p{p - Q) . (29) 
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Of course we also have Vp{z) = in such a way that the generic primary VpVp has 

dimensions (Ap,Ap). 

We have already seen that at criticality the requirement of mutual locality of the fields with 
the energy-momentum tensor is implemented by (|20p . When moving to the off-critical action 
(HOD, the energy-momentum tensor acquires a non-zero trace proportional to the field responsible 
for the breaking of scale invariance: Q{x) oc t£{x). We now examine, first for c = 1 and then for 
c < 1, some implications of the choice of e and of the requirement of locality with respect to e. 


3.2 c=l 

Central charge equal 1 corresponds to the case Q = 0, namely to the usual Gaussian theory 
with OPE 

^Pl ■ ^P2 ~ [^1-|-P2] 5 (30) 

where, with respect to (HZD, we omit for brevity coordinate dependence and structure constants; 
the brackets indicate that, besides the primary field, the r.h.s. contains its descendants. Using 
(fT8l) . (I29D and (l30D we obtain that two fields Vp-^Vp^ and Vp^Vp^ are mutnally local if 

2(piP2-P 1 P 2 ) e Z. (31) 


A generic choice for a real and scalar energy density field in this theory is e = I4I4 -|- V-bV^h oc 
cos2bip with Aj = 6^; a field VpVp is local with respect to e if ±26(p — p) is an integer, i.e. if 

P-P = ^, meZ. (32) 

Taking m = 1, we can build the complex (Dirac) fermion 


T = (V^,V') 


V±,bV_^,b,V±_bV_^_b] , 

46 “'"2 46 “'"2 46 2 46 2 / 


(33) 


with spin equal 1/2 for and —1/2 for pj] the decomposition T = Ti -|- i '^2 defines two 

real (Majorana) fermions Tj = For 6^ = 1/2 we have dip = dPp = 0, which are free 

fermionic equations of motion; hence, for 6^ = 1/2, the theory (1231) admits a representation in 
terms of free fermions. For 6^ 7^ 1/2, instead, the fermions are coupled by the four-fermion term, 
which can be shown to be truly marginal, so that the bosonic action (I23p with energy density 
held cos 2hip fermionizes into [mHZKIl] 


An = drx 




{'ipidipi + + g{b^)ipi'ipi'4’2ip2 


i=l,2 


(34) 


with g{l/2) = 0; cos26(/? fermionizes into the mass term 'ipi'ipi -|- V'2^2- 

Since the real fermionic helds satisfy ipf = 'ipf = 0, so that ^pi'4>ii^2ip2 oc (X)jthe 
action (p^ is invariant under the 0(2) rotations of the vector (V’i,V'2); these become U{1) 
transformations for the complex fermion, and the integer m in ()32p plays the role of U{1) 
charge. Then the components of the order held associated to the 0(1) ~ 0(2) symmetry 
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are selected picking up the scalar fields with m = ±1, i.e. a± = V"±i/4feKpi/4b with = 
1/166^. If, by construction, cos26(/9 is the most relevant 0(2)-invariant field, generic fields 
entering the most general action with 0(2) symmetry must have m = 0 and must be local 
with respect to the order field components a±; ([3T]) then identifies these fields with with 

k integer and conformal dimension k'^b'^. We then see that, ignoring the identity {k = 0), 
all these fields are irrelevant for 6^ > 1; cos2bip is the only invariant relevant field for 1/4 < 
b^ < I and leads to an off-critical action (known as sine-Gordon model) which, as we will see 
later, is purely massive for both signs of the conjugated coupling. This leaves no room for a 
phase with spontaneously broken 0(2) symmetry, which requires a relevant invariant field and 
massless (Goldstone) excitations associated to a continuum of degenerate vacua. Indeed, the 
two-dimensional XY model only exhibits a transition between a high-temperature phase with 
exponential decay of correlations and a low-temperature phase with algebraic decay, while the 
order parameter vanishes at all temperatures. This is known as Berezinsky-Kosterlitz-Thouless 
(BKT) transition [TO] and is described by the Gaussian model with b‘^{T) a (non-universal) 
decreasing function of the temperature; 6^ = 1 at the BKT transition temperature Tbkt, so 
that the low-temperature regime with algebraically decaying correlations is naturally explained 
as the regime in which all 0(2)-invariant fields are irrelevant and the theory is conformal at 
large distances. Since = 1/16 at b^ = 1, spin-spin correlations decay as at Tbkt- 

The XY model provides the basic illustration of the general fact that continuous symmetries do 
not break spontaneously in two dimensions |2ni[2ll[22|. 

It is known from exact lattice results |23] that the two-dimensional Ising model corresponds 
to the theory of a free real fermion. Hence, the action (1341) with 6^ = 1/2 describes the scaling 
limit of two decoupled Ising models, each with central charge c = 1/2, energy density e* = 
and order field cTj {i = 1,2). The conformal dimensions Ag = 1/2 and Ao-j_ = 1/8 at this 
decoupling point coincide with the dimensions A^. and A^^^^ = 2A^^, respectively, which follow 
from exact results for the lattice Ising model [241125j . 

For later use we are also interested in the chiral (i.e. depending on z or z only) fields which 
satisfy (f32]l and have lowest charge m = ±1; they are rj± = Y±\i 2 h aiid = ^±i/2b) with 
Ajjj_ = Afjj_ = 1/46^. For b generic the spin of these fields is neither integer nor half-integer, and 
for this reason we call them “parafermionic” fields isaEZ). 

3.3 c<l 


In the case of central charge smaller than 1, 

which corresponds to Q real, it is useful to param- 

eterize Q as 



Q = /3-i-/3, 

— t/{t + 1), 

(35) 

and to introduce the notation 





(36) 
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where = ^[(1 — ^ — (1 — i^)/3]. Then the central charge ([27)1 can be rewritten as 


c = 1 — 


and 


t{t + l) ’ 

gives for the single non-vanishing conformal dimension of the fields (|36[) the result 

[(t + l)/i - - 1 


(37) 


= 


4t(t + 1) 


(38) 


It is known from [12] that when the indices and v take positive integer values m and n the 
dimensions ([55|) are those of the degenerate primary fields associated to the algebra (ITTIl for 
c < 1, and that the differential equations satisfied by multi-point correlators with degenerate 
fields imply the degenerate-non-degenerate and degenerate-degenerate OPE’s 


m—1 n—1 

^m,n * ^ ^ ^ ^ ^ [^//—m+l+2fc,t/—n+1+2/] ? (^9) 

k =0 1=0 

min(mi,m2 )—1 min(ni ,712) —1 

* ^m 2 ,n 2 — E E m2|-|-l-|-2fc,|ni—n2|+l+2/] ! ( 4 ^ 0 ) 

k =0 1=0 


with similar relations for the fields 

The energy density field e is consistently defined if the OPE e ■ £ does not produce fields 
(other than the identity) with scaling dimension smaller than and, within the present 

knowledge of theories with c < 1, this seems to require that e is a degenerate field. We further 
restrict our analysis to theories which reduce to the Ising model for c = 1/2, i.e. have = 1/2 
for t = 3. Then (l38l) and (170l) leave us with $2,1^2,! and 3 as possible choices for s. Eor 

the two choices, the OPE’s ([39]l and (I70l) can be used to select the fields mutually local with e 
and, among these, the parafermion r]{z) and the order field which has to satisfy 


£ ■ a = a + ■ ■ ■ 


(41) 


on symmetry grounds. 

Eor £ = <^2,1^2,1 this analysis |25| yields rj = <^1^3 and a = ^i/ 2 ,o^i/ 2 ,o- In this case the OPE 
rj ■ a produces a field /r with the same conformal dimensions of cj; the reason we do not identify 
with a is that we expect r] to be charged with respect to the symmetry G (to be identified), 
as we saw for c = 1, so that a and cannot transform in the same way. The field is called 
“disorder” field and its presence is the signature of an order-disorder “duality”. Eor e = $1,3^1,3 
one finds r/ = $2,1) while (|4ip does not help for the identification of a. 

The two choices of £ define critical lines parameterized by t in the space of c < 1 conformal 
field theory. Before discussing them in the next section, we mention that c < 1 conformal field 
theory also allows for discrete series of critical points with special properties. Indeed it can 
be seen that for rational values of the parameter /3^ in ()35l) the OPE (I40p closes on a finite 
number of degenerate primaries; one can then build consistent conformal theories containing a 
finite number of primaries and their descendants which are known as “minimal models” [12] . 


14 




Figure 3: Pictorial illustration of equation (14211 in (1 + l)-dimensional space-time. 


In particular, for t = 3,4,..the conformal dimensions Am,n of the primary fields are given 
by (I38p with m = 1,... ,t — 1, n = 1,... ,t, and are all non-negative [29]; this set of minimal 
models describes the critical {t = 3) and multicritical {t > 3) points associated to Z 2 symmetrj|§ 
isniEi]. The model with t = 5, restricted on a smaller set of primaries [321133] . corresponds to 
the critical point of the three-state Potts model |34] . as we will also see in the next sectior0. 


4 Scale invariant scattering and critical exponents 

4.1 Scattering and statistics 

We have seen in section [2] that the canonical critical exponents are determined by the scaling 
dimensions = 2A^ and Ao- = 2Ao- of the energy density and order fields. We also introduced 
the g'-state Potts and n-vector models and saw that their continuations to real values of q and n 
are related to percolation and non-intersecting walks, respectively. Since the Ising model, which 
corresponds to g = 2 and n = 1 , has central charge c = 1 / 2 , these continuations should produce 
lines of critical points contained, entirely or in part, within the space of theories with c < 1 
discussed in the previous section. In order to identify A^ and A^- in the two cases we need to 
relate these conformal field theories with the symmetries Sg and 0 (n) which characterize the two 
universality classes. Following |28|, we will do this using the particle description introduced in 
section 1231 The results we will obtain for the exponents coincide with those originally determined 
within the lattice Coulomb gas framework |3H] and cast in the framework of conformal held 
theory in m- 

The elementary excitations of a scale invariant relativistic (l-l-l)-dimensional theory are 
massless particles with energy and momentum related a^ = e > 0 (right movers) or p = —e < 0 
(left movers). These particles are created by the chiralj helds rj, with Arj = 0, if right movers, 

^The scaling dimensions of primaries in the minimal model with t = 3 (Ising) are Aip = A 2,3 = 0, Ai ,2 = 
A 2,2 = 1/16 and Ai _3 = A 2 ,i = 1/2, corresponding to the identity, a and e, respectively. 

“^We also mention that the precise status of degenerate fields is the object of current studies in the so-called 
“logarithmic” cases m, where the coincidence of the dimensions of two fields may lead to logarithmic corrections 
to the algebraic behavior (0 of correlators (see [36] for a recent review and a list of references). Normally, the 
collision of scaling dimensions and the associated logarithmic corrections can be observed at special points on 
critical lines with ordinary power law correlations (see e.g. |37p. 

® Recalling da we have (0|4>(2:)|p) proportional to for a right mover and to for a left mover. 
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Figure 4: Amplitudes 5i, S 2 and ^3 of the 0(n)-invariant theory. 


and f], with Af^ = 0, if left movers; one also has Sr^ = —Sf^ = = Af^ for the spin of the fields. 

Due to the absence of dimensional parameters, the dimensionless elastic scattering amplitude, 
that we denote by 5, of a right-mover with a left-mover cannot depend on the variable (jlljl . 
which is the only relativistic invariant in the process and is dimensionful. The fact that S does 
not depend on the particle energy and, by unitarity, is a constant phase means that the two 
particles have no dynamical interaction. The scattering process, however, involves exchanging 
the position of the two particles on the line, so that scattering entails in general a statistical 
factor; S will be 1 for bosons (A^ integer), —1 for fermions (A^ half-integer), and a more general 
phase when the helds rj and f] are parafermionic. Since in absence of dynamical interaction the 
passage from the initial to the hnal state can also be realized by vr-rotations (see Fig. [3]) ruled 
by the Euclidean spin, the right-left scattering (statistical) phase can be written as [28] 

5 _ g-iTr{Sr,-Sfj) _ ^-2iTTAr, _ /^ 2 ) 


In general the particles and the associated chiral helds form multiplets transforming under a 
representation of the symmetry of the theory. We denote by r]a and fia the components of such 
multiplets, and encode the general right-left scattering in the symbolic relation 


VaOVb = Y^ Sat VdOVc, 

c,d 

where S^t are scattering amplitudes (Fig.[T|); then (|4^ is obtained upon diagonalization of 
Unitarity yields 


(43) 


cd 


e,/ 

while crossing symmetry ()12p can be written as 




ocd _ 

^ab — 


r>bc 

"^da 


(44) 


(45) 


recalling real analyticity. We can now proceed to the analysis of specihc cases. 


4.2 0{n) symmetry and non-intersecting walks 

In this case the symmetry is naturally represented by vector multiplets rja, fja, a = 1,... ,n, (jSj) 
takes the 0(n)-covariant form (Fig. |4|) 

n 

r]aOVb = Sab 5i ^ fyc O rye + S 2 % o -f S 3 fjaorib, (46) 

C=1 
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gives 5i = ^3 = /Oie**^, 52 = = ^2, with pi non-negative and p 2 reaO, and is then 

rewritten as 


Pi + /’i = 11 (47) 

P1/32COS99 = 0, (48) 

np\ + 2pip2 cos p + 2p\ cos 2(/? = 0 . (49) 

There are three ways of satisfying (1481) . The first possibility is that cosp = 0, so that (I49p fixes 
re = 2, a case to which we will come back in a moment. The second possibility is that pi = 0, 


but then we are left with 52 = ±1 as the only non-vanishing amplitude, i.e. with re decoupled 
free bosons or fermions. Hence, the only non-trivial case with a continuous re-dependence is 

P 2 = 0, pi = 1, re = —2 cos 2p S (—2,2). (50) 


Thinking of particle trajectories as walks of re different colors on the Euclidean plane, the result 
52 = 0 amounts to a no-crossing condition for the walks (Fig. H]), as Hrst observed in [39] in the 
study of the off-critical case. 

It follows from ([46]l and (l50l) that 


n n 


^ Pa 0 Pa = 5 ^ Pa 0 Pa , 
a=l a=l 

(51) 

5 = re5i + 52 + 53 = , 

(52) 


and (|42l) then relates p to up to the 2TTk ambiguity involved in the comparison of phases. 
The relations (I5UI) and (I52p consistently give re = 1 and the fermionic (Ising) value 5 = — 1 
when p = 2'Kkl‘i^ k = 1,2 (mod 3). Looking for matching with one of the two c < 1 critical 
lines identified in the previous section, we take that with e = 4>i^34>i^3 and p = $2,1, since the 
order-disorder duality exhibited by the other is typical of the Potts model mm- Then we use 
Aj = Ai ^3 = 1 at c = 1 and Aj^ = 1/4 for Ag = 1 (i.e. for 6^ = 1) on the 0(2) line of the 
previous section to fix = —^(A^ -|- 1/2) (mod27r); substituting A^ = A2,i and plugging into 
(f50P we finally obtain 


TT 


re = 2 cos — 
t 


(53) 


where t parameterizes the central charge through (|37p . We see that the re = 2 endpoint of the 
solution (I50p corresponds to the BKT transition point 6^ = 1 on the 0(2) line with c = 1 of the 
previous section. Hence we conclude that (l50P gives A^- = 1/16 not only at c = 1/2 but also 
at c = 1; using these conditions to determine the indices (assumed t-independent) in ()38p gives 
Act = A 1/2,0 S’® the only positive solution within the range c G (0,1). 

Coming back to the solution of (|47II49I) with re = 2, it can be parameterized as 


Pi = sin a, p2 = cosa, (p = —7r/2, 

®The phase p should not be confused with the field in (1231) . 


(54) 
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Figure 5: Amplitudes Sq, 5i, S 2 and ^3 of the S'^-invariant theory. 


where we now allow for pi negative. The phase in (|5ip becomes 


5 = 25i + 52 + 53 = e-*“ 


(55) 


and has to be equated to ()42p with which, as seen in the previous section, takes the value 
1/46^ along the 0(2)-invariant critical line parameterized by 5^; this gives 


TT 


(56) 



a relation that we will exploit later on. 

4.3 Permutational symmetry and cluster boundaries 

In view of the cluster representation ([7|) of the Potts partition function, we associate in this 
case the particle trajectories to the boundaries separating different clusters, denote by fjap 
(a, /? = 1 ,..., g; a 7^ /3) the chiral helds corresponding to the trajectories separating a cluster of 
color a from one of color (3, and use the permutational symmetry of the colors to write (1431) in 




(57) 


Crossing symmetry enables us to write 


(58) 



with po and p3 real and p non-negative, so that the unitarity equations take the form 


{q - 2 >)pI + p^ = 1 , 

{q - ^)Po + ‘^PoP cos = 0 , 
{q — 2)p^ + p^ = l, 

{q - 2>)p^ + 2pp3 cos = 0 . 


^This color structure for the scattering problem was first observed in [40] within the study of the model below 
Tc, with the particles associated to kinks interpolating between degenerate ferromagnetic vacua. 
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Symmetry 

c 

Ae 

A^ 

A. 

0 ( 2 ) 

0 (n), n = 2 cos j 

Sq, ^ = 

1 

1 ® 

1 ® 

t(t+i) 

Ai,3 

A2,1 

1 

iP' 

A2,1 

Ai,3 

1 

TeP 

^ 1 / 2,0 

^ 1 / 2,0 


Table 1: Central charge and conformal dimensions along the three critical lines discussed in the 
text; is given by (f38|) . The dimensions and Ao- of the energy density and order helds 
determine the canonical critical exponents. 


The solution 

/>o = -l, p = - g, 2cos(p = = q - 3 ^ (59) 

is the on^ which identifies a S'q-invariant critical line with q £ (0,4), in agreement with the fact, 
known exactly from lattice studies m, that the Potts phase transition becomes of the first order 
for (? > 4, thus preventing a continuum limit. At q = 2 the only physical amplitude is ^3 = — 1, 
as needed for the Ising model, and the solution (I59p must be associated to the remaining theory 
with c < 1 of the previous section, that with e = 4 > 2 , 1 ^ 2,1 and rj = c = 1 should correspond 
to the maximal value of g, which we saw is 4. It follows from (1571) and (|59p that 

^ ^ ® P 70 — S ^ ^ 'IJa'y ® V'ya ; (60) 

7 7 

S = S3 + {q- 2)52 = . (61) 

Using (|i^ and the requirement that A^ = Ai ^3 equals 1/2 and 1 for 5 = 2 and 4, respectively, 
one finds (/? = 7r(l + ^Ai^ 3 ) (mod27r) and, upon substitution in (l5^ . 

^ 7r(t — 1) , , 

Table 1 summarizes the results for the three critical lines we have discussed in this section. 

5 Away from criticality 

5.1 Integrability 

The fact that two-dimensional scale invariant theories yield exact results (in particular the 
critical exponents) has to be traced back to the infinite-dimensional conformal symmetry codified 
by the algebra m, which in turn implies the existence of infinitely many conservation laws. A 
conservation equation in two dimensions takes the form 

BTs+i = dQs-i , (63) 

®We use the condition pa = — 1 at g = 2 (Ising) to fix the sign of ps. 
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where the subscripts specify the Euclidean spin of the fields. At criticality, the descendants 
of the identity with conformal dimensions (A, A) = (s + 1,0) yield infinitely many local fields 
Ts+i satisfying (|63p with zero on the r.h.s; hence, that of energy and momentum (s = 1, 
Eq. (llOp l is only one among infinitely many conservation equations. When we move away from 
criticality and consider a theory with action m, dTs-\-i no longer vanishes, but for s / 1 
has no need to become a total derivative with respect to z. In other words, non-trivial (i.e. 
other than energy and momentum) conserved quantities are normally lost away from criticality. 
Instead of (1631) we will have in general the expansion [l2] dTg^i = '^^=1 ^ where r with 

positive dimensions (1 — A^, 1 — A^) is the coupling in (llOp . and are fields with dimensions 
(s + 1 — n(l — Ag), 1 — n(l — Ag)). We consider theories with a spectrum of dimensions discret^ 
and bounded from below. Then N not only has to be finite, to avoid fields with arbitrarily 
large negative dimensions, but generically has to be 1, since only = Ag can be identified 
with a field certainly present in the theory, i.e. the descendant of £ with dimensions (A^ + s, Ag), 
without special conditions on A^. Hence, in the generic case we are left with 

dTsJ^i = TAs, (64) 

which is a conservation equation if Ag is d of another field; of course we exclude the case in 
which Tg^i is itself d of another field. Let us call the number of linearly independent non¬ 

derivative fields Tg^i, i.e. the dimension of the space of level s -|- 1 descendants of the identity 
which are not d of other fields, and the dimension of the space of level s non-derivative 
descendants of e; then (I64h is a conservation law if 

Das < Dt,+i ■ (65) 

Such a counting argument |42] yields a sufficient condition for the existence of off-critical con¬ 
served quantities through the comparison of dimensionalities of subspaces in the space of fields; 
the structure of these subspaces has been described in section 13.11 Eor s = 1 we have Dai = 0 
(the only level 1 field is L_ie = de) and Dt2 = 1 {L-2I = T), so that (1651) is always satisfied 
and conservation of energy and momentum recovered. Eor s > 1, as expected, (|65p is violated 
in general. For example, for s = 3 one has Da^ = = 1, since L-^e and the 

only fields allowecj^. However, if e is a primary degenerate (in the sense defined in section 13.11) 
at level 3, Da^ reduces to 0 and (1651) is satisfied. Hence we see that, remarkably, if the field 
which breaks conformal invariance is degenerate, the reduced dimensionality of the descendant 
subspaces leaves room for ()65p to be fulfilled with s > 1, i.e. for the existence of non-trivial con¬ 
served quantities. The systematic analysis performed in [42] shows that the relevant degenerate 
fields </>i, 3 , 01,2 and 02,1 allow for non-trivial conserved quantities when one of them is added to 

®While the assumption of discrete spectrum of dimensions covers the main cases and is certainly sufficient for 
our purposes, it is interesting to mention that lattice models pointing to a continuous spectrum of dimensions 
have been investigated |43| . 

^°The derivative character of L-nI follows from the fact that (|20)) with <1? = 7 reduces to a Taylor expansion 
for T{z). 
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Figure 6: Factorization of scattering amplitudes in integrable field theories. 


a conformal field theory with c < 1; here we are introducing the notation 

4>im,u{x) = ( 66 ) 

for scalar primaries obtained from the fields (j36p . 

While the sufficient condition ()65p turns out to yield a finite number of conservation laws, 
it is expected that these are a subset of an infinite series which leads to the exact solvability (or 
integrability) of the off-critical theory. Ultimately, this expectation is confirmed by the existence 
of exact S'-matrix solutions, as we are going to see. Anticipating this step and recalling Table [H 
we observe that the above results of the counting argument mean that the two-dimensional q- 
state Potts {q G (0,4)) and n-vector (n G (—2,2)) models are exactly solvable in the scaling limit 
away from criticality. The same conclusion follows for the Ising model at T = Tc in presence of a 
magnetic field (the magnetic field amounts to the addition of H ai to the lattice Hamiltonian 
(|3|)). The scaling limit at Tc is obtained adding to the conformal theory with central charge 
c = 1/2 the field a = ^ 1 / 2,0 = 0i,2; the last equality follows from (p8]l . which for t = 3 gives 
^ 1 / 2,0 = '^ 1,2 = 1/16. This is one of many examples of models which are not solved on the 
lattice but become solvable in the scaling limit. 

If a field theory possesses inhnitely many conserved quantities its 5-matrix drastically sim- 
plihes. First of all the infinitely many conservation equations force the final state of a scattering 
process to be kinematically identical (same number of particles, same masses, same momenta) 
to the initial one [complete elasticity)] only reshuffling of charges is allowed. The second sim¬ 
plification can be understood through the following argument [H]. A single-particle state \p) 
is eigenstate of a spin s conserved quantity with eigenvalue proportional to [p + e)^. If the 
particle is located around xi = 0 at t = 0 it can be described by the wave packet f dkf[k)e^^^^, 
with f[k) peaked around p. Stationary phase approximation then shows that the action on 
the wave packet of the translation operator moves it to xi = a, no matter the value of 

p. If, on the other hand, we repeat the operation replacing the momentum operator P with 
the spatial component of a spin s conserved quantity, this produces a factor and the sta¬ 
tionary phase approximation now yields a displacement proportional to p^~^. This means that 
non-trivial conserved quantities (s > 1) act on particles as generators of momentum-dependent 
spatial displacements. This fact can be used in a multiparticle scattering to make arbitrarily 
distant in space-time the two-body subprocesses, as well as to change their chronological order 
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Figure 7: The cut s-plane is mapped onto the cut-free 0-plane through the relation ()69l) . 


(Fig. [6|). Since the conserved quantities commute with time evolution, we conclude that n- 
particle scattering amplitudes factorize into the product of n{n — l)/2 two-particle amplitudes, 
and that 

5i25i3523 = 5235i35i2 , (67) 

where Sij is a short notation for the amplitude of the particles with momenta pi and pj. Let us 
stress that Fig. El and then factorization, only makes sense in 1-|-1 dimensions, since two particles 
with different momenta necessarily meet on a line. As we will see below, when added to unitarity 
and crossing symmetry, complete elasticity and factorization allow the exact determination of 
the S'-matrix. 

Factorization of scattering amplitudes in integrable field theories reduces the determination 
of the S'-matrix to that of the two-particle amplitudes of Fig. □ [15]. Complete elasticity of the 
scattering implies that these amplitudes exhibit in the s-plane only the elastic unitarity cut and 
its crossing-symmetric image (Fig. [7]). It is convenient to introduce the rapidity parameterization 

(e,p) = (mcosh0,msinh0) (68) 

for the momentum and energy of a particle of mass m; a Lorentz transformation corresponds to 
a shift of the rapidity, so that relativistic invariant quantities depend on rapidity differences. In 
particular, the square of the center of mass energy (ttH) becomes 

s = 4m^ cosh^ (69) 

for two particles of equal mass and rapidities 0 i, 62 ( 0 i 2 = 0i — 62 ), and we will use the notation 
~ amplitudes. The fact that different values of 9 correspond to the same 

value of s can be exploited to map different sheets of the s-plane onto different regions of the 9- 
plane, the physical sheet corresponding to the physical strip Im0 G (0, vr) (Fig.[7|); the amplitudes 
S^f{9) are then “meromorphic” functions, i.e. allow for poles as the only singularities in the 
complex 0-plane. Due to complete elasticity, the unitarity condition (fTHI) holds at any energy, 
and in the rapidity variables becomes 

Y,s:l{9)Sff{-9) = dac 6 bd; (70) 

ej 
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Figure 8: A bound state pole of the scattering amplitude. 


the crossing relation (1121) translates into 


(71) 


5.2 Purely trasmissive scattering 

When looking for solutions of the factorized scattering problem the simplest possibility is to 
consider purely transmissive scattering, i.e. c = o and d = 6 in Fig. [TJ For this case we can set 
= Sab{d) and (ITOl) . d7T|) reduce to 


Sab{0)Sab{-0) = 1 , (72) 

Sab{e) = S-,,{i7r-e). (73) 

The solutions to (j72l) in the space of functions meromorphic in 0 and real analytic and polyno- 
mially bound in s can be written in the form n«/a('9)> with 


sinh[(d + i'Ka)/2] 
sinh[(0 — i'Ka)/2] ’ 


a G (0, 2). 


(74) 


A pole 

Sah{0) ^ ^ - i'^ab (75) 

^ - Kb 

in the physical strip corresponds to a bound state with mass rric = 2mcos(uJ(j,/2) (Fig. [8|). Up 
to the factor of i, such a direct channel pole has a positive residue; through ()73p it induces a 
negative residue pole in the crossed channel. 

Trivial examples of purely transmissive scattering are provided by free theories. It follows 
from the results of section [3] that these include the relevant case of the scaling Ising model which, 
in absence of magnetic field, corresponds to the theory of a free massive real fermion, and then 
to S{9) = —1. This amplitude describes both scaling regimes above and below Tc, since the 
difference in the action is only a change of sign of the mass term r'0^. Below T^, however, 
there are two degenerate vacua and the free fermions corresponds to the kinks which interpolate 
between them, which in (1 + 1) dimensions are stable for topological reasons and always provide 
the elementary excitations in presence of degenerate vacua. If denotes a kink going from 
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Figure 9: A bootstrap equation. 


the vacuum a to the vacuum /3, in the Ising case we only have _and K _and multi-kink 

states can only have the form • • • _ K __• • •, with the sequence fixed by the choice of the 

first vacuum. By symmetry there is a single two-kink amplitude, which can be depicted as the 
last one of Fig. [5] and is equal to —1. The coincidence of the amplitudes above and below Tc is 
the manifestation, within the particle description, of the duality of the Ising model well known 
from the lattice formulation |46j . 

To illustrate the interacting case we consider a theory invariant under a charge conjugation 
symmetry and whose elementary excitations are two conjugated particles that we denote by A 
and A] they have charge q and —g, respectively. The simplest non-constant solution of (1721) and 
(1731) reads [17] 

SAA{e) = f2/3{0) , S^AiO) = -fl/3i0) • (76) 

The overall sign is chosen so that the pole at 6 = 2i7r/3 of Saa is positive, and then corresponds 
to a bound state with the same mass m which does not introduce new particle species. Identifying 
this particle with A would imply g = 0; it then corresponds to A, and we have 2 q = —q, i.e. 
identification of charge modulo 3q. Hence, the amplitudes (|76p yield the exact solution of the 
simplest theory possessing symmetry {S3 taking into account also charge conjugation). The 
freedom to separately shift particle trajectories in integrable theories also yields the “bootstrap” 
equation (Fig. [9|) 

'S'aa(^) = SAA{d - in/3)SAAid + W3) (77) 

in connection with the vertex A A —>• A; it is satisfied by ()76p . 

On universality grounds it is natural to identify (|76l) to the scattering solution of the scaling 
three-state Potts model, whose integrability we argued above. The identification is correct 
[481 ll9] and in this as in the other cases the ultimate check is provided by a technique known 
as thermodynamic Bethe ansatz |l9]. In this approach one studies the thermodynamics of the 
gas of relativistic particles confined on a circle, and uses the scattering amplitudes to specify 
the interactions in the low-density limit. When the temperature of the gas becomes much larger 
than the mass of the particles the theory approaches the scale invariant limit and it is possible 
to determine the central charge starting from the scattering amplitudes. 

The above solution in terms of A and A corresponds to the scaling three-state Potts model 
above Tc. Below Tc there are three degenerate vacua and the associated kinks, together with 
the symmetry, allow for the inequivalent amplitudes Si, S2 and S3 corresponding to the last 
three in Fig. [5l The particle-kink correspondences A -H- Ka^a+i, A -H- with indices 
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taking integer values identified modulo 3, allow the identifications Saa = *S'i, S^a — ^2, and 
then the description of both regimes above and below in terms of the same amplitudes 
(53 = = 0). This is expected, since the lattice g-state Potts model exhibits high-low 

temperature duality mm- 

As another relevant example of purely transmissive theory we mention the Ising model at 
Tc with a magnetic field, whose integrability we also discussed. The presence of the magnetic 
field leaves no internal symmetry, so that the mass spectrum has to be non-degenerate and 
the scattering transmissive. A detailed analysis of non-trivial conservation laws at bound state 
vertices was exploited in [32] to determine the full particle spectrum (which consists of eight 
particles with different masses) and the scattering amplitudes. The reader is referred to [50] for 
a review of this solution and of its implications. 

5.3 In presence of backscattering 

The case in which the scattering is factorized but not purely transmissive can be illustrated 
through the very relevant example of the sine-Gordon model. This corresponds to the action 
(|10l) with the conformal invariant part App given by the c = 1 free massless bosonic action 
(ESD, and e oc cos2bip. It follows from the discussion of section that the model also admits 
the fermionic representation with ()34p as the conformal part and e oc (V’lV'i + '</’ 2 V’ 2 )- We also 
discussed how the complex fermion (I33p and its conjugate T* carry opposite units of the charge 
associated to the U{1) symmetry; these helds create the elementary particle excitations of the 
theory, that we denote by A and A, respectively. In the bosonic language they correspond to 
the soliton and antisoliton which interpolate, in opposite directions, between adjacent minima 
of the periodic potential cos2bip-, the 17(1) charge then corresponds to a topologic charge. 

The integrability of sine-Gordon model does not follow from the counting argument, which 
does not apply to c = 1, but can be shown using more traditional methods of Lagrangian field 
theory (|35| and references therein). Hence the scattering solution requires the determination 
of the amplitudes S^f{9) with indices taking values A and A. Charge conservation and charge 


conjugation symmetry allow for the three independent amplitudes 

Sii{9)^So{e), Sij{9)^ST{9), sii{9)^SR{9), (78) 

which satisfy the crossing equations 

So{9) = SriiTT - 9 ), 5^j(0) = 5jj(z7r - 9) , (79) 

and the unitarity equations 

STi9)ST{-9) + Sr{9)Sr{-9) = 1, (80) 

STi9)SR{-9) + Sr{9)St{-9) = 0, (81) 

So{9)So{-9) = l. (82) 
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The presence of the reflection amplitude Sr makes non-trivial the factorization equations (I67p 
since, for fixed initial and final states, one has to sum over the different particle indices possible 
for the internal triangles in Fig. [6l In particular, the process AAA —)■ AAA gives 

ST{9)So{e + 9')SR{e') + Sr{9)Sr{9 + 9')St{9') = SR{9')STi9 + 9')So{9) . (83) 


For the theories in which the scattering is not purely transmissive the amplitudes cannot be 
written as products of basic blocks like (1741) . and their determination requires a non-trivial 
model dependent study of the functional equations. For sine-Gordon this can be found in [45j 
and leads to the minimal (i.e with the minimal number of poles) solution to (I79II83P 
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the amplitude Sq is initially found in the form of an infinite product of gamma functions, more 
suitable for dealing with the functional equations, and then exponentiated using the formula 


F( 2 ;) = exp 
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( 86 ) 


The parameter ^ allowed by the solution of the functional equations for the amplitudes must be 
related to the parameter b which in the sine-Gordon action provides the coordinate along the 
0(2)-invariant critical line of section 13.21 The amplitude Sr vanishes and Sr and Sq become 
— 1 when ^ = vr, which then corresponds to the free fermion point 6^ = 1/2. The limit 0 —>• oo 
(i.e. m/v^ —)■ 0) of Sr, which is e 2 V should coincide with the massless amplitude ([5^ . 
Using (I56|) and the condition at ^ = vr one obtains the relation 




7r6^ 


(87) 


originally obtained from semiclassical and perturbative studies ( [l5] and references therein). The 
requirement that the energy density cos26(/3 is a relevant field implies 6^ < 1, so that ^ > 0. 
For ^ < TT the amplitudes Sr and Sr exhibit simple poles in the physical strip corresponding to 
soliton-antisoliton bound states Bn. The poles are located at 9 = i{'K — nQ and determine the 
masses of the particles Bn as 


rrin 


o ■ ni 
2m sm — , 
2 ’ 


1 < n < 


TT 

? 


( 88 ) 


where [x] denotes the integer part of x. The lightest bound state Bi is the particle created 
by the field ip; it is no longer in the spectrum of the theory for ^ > vr, which is a regime of 
repulsive soliton-antisoliton interaction. The scattering amplitudes involving the particles Bn 
can be obtained from those for A and A exploiting bootstrap equations similar to that depicted 
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in Fig. [9l the poles of these new amplitudes show, in particular, that Bn can be seen as the 
bound state on n particles Bi (see [IS])- 

The fermionic form of the model corresponds to the scaling limit of two Ising models coupled 
through the product of their energy densities (the four-fermion term in dMj))- This is known as 
the Ashkin-Teller model and is defined on the lattice by the Hamiltonian 

'Hat = — -I- a2^i(T2,j\ + , (89) 

(b> 

where ai^i and < 72,1 are the two Ising spins at site i. This four-state model possesses a line of 
second order phase transition in the space (J, J 4 ) which is described by the line of conformal held 
theories with c = 1 , and was originally cast in the bosonic/fermionic held theoretical framework 
in [SUES]. For J = J 4 the Ising variables a 2 ^i and <ti play a symmetric role in (f89]l . 
corresponding to the four-state Potts model subspace. On the other hand, it follows from Tabled] 
that the critical four-state Potts model has c = 1 and = 6 ^ = 1/4. As a consequence, the sine- 
Gordon scattering amplitudes with ^ = tt/S solve the scaling four-state Potts model; notice that 
Sr vanishes at this point, so that we have purely transmissive scattering as in the three-state 
case, and one can check that the remaining amplitudes reduce to products of the blocks (1741) 
(the same is true for all integer values of tt/^). At the particle level the symmetry enhancement 
is signaled by the fact that for = tt/S the mass of the particle Bi becomes degenerate with 
that of A and A. More generally, the scattering description of the Ashkin-Teller model is given 
in [531 [M] and we will come back on this in section 01 

5.4 Antiferromagnets 

An antiferromagnetic interaction corresponds to taking J < 0 in a lattice Hamiltonian such 
as m, so that configurations in which nearest neighboring spins take ‘opposite’ values are 
energetically favored. Since the extent to which spins can differ from their neighbors depends on 
the number of neighbors, phase transitions in antiferromagnets are crucially dependent on the 
lattice structure and, in this sense, non-universal. Still, if an antiferromagnet on a given lattice 
undergoes a second order transition, this will be described by a held theory in the scaling limit. 

Considering for simplicity the square lattice, the Ising case is particularly simple, since it 
can be mapped back to the ferromagnetic case introducing the staggered spin variable Sj = 
(_l)n+* 2 ^^, where i = ( 11 , 12 ) is the site label. The three-state Potts model provides instead a 
non-trivial illustration of the combined effect of internal and lattice symmetries. Now there are 
inhnitely many conhgurations in which each spin has a value different from that of its nearest 
neighbors, so that the model possesses inhnite ground state degeneracy. These ground states 
admit an alternative representation in terms of arrows on the edges of the square lattice. After 
transferring the labels (colors) Si = 1,2,3 from sites to faces, an observer in a face with color j 
looking across an edge to an adjacent face with color j ± 1 (mod 3) puts an arrow on this edge 
pointing to his left/right. The arrow configurations obtained in this way satisfy at each site the 
rule “two arrows in, two arrows out” of the six-vertex model, an exactly solvable model [551I1I] 
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which is critical and is described in the continuum by c = 1 (Gaussian) conformal field theory. 
Hence the three-state Potts antiferromagnet on the square lattice is critical^ at T = 0; this in 
turn means that its scaling limit for T > 0 is described by the sine-Gordon model with energy 
density e oc cos 2bip, for some specific value of b. 

In order to further understand the nature of the correspondence with the Gaussian model, we 
consider the complex lattice spin variable aj = staggered one Sj = 

where Sj = 1,2,3. Gontrary to the energy density, Sj and aj carry a charge associated to 
cyclic color permutations; Sj is also odd under exchange of even and odd sublattices, namely 
the two sublattices identihed by the parity of ji + j 2 - Hence in the scaling limit these variables 
correspond to charged scalar fields of the sine-Gordon model, in which charge is counted in units 
of the integer m of ([32]) . Using the notation Um{x) = Pm/ 46 ( 2 )U_m/ 46 ( 2 ) for the charged scalar 
fields, the natural identifications are S = C/i, o" = U -2 = U 2 , consistent with the identification of 
the Z 3 charge with m (mod 3), and of the sublattice parity with (—1)”* [57]. For the conformal 
dimensions we have As = 1/166^, A^- = 1/46^ and A^ = 6^; they match the lattice results (see 
|58j and references therein) for 6^ = 3/4. The amplitudes (|84l) with this value of 6^ yield the 
scattering solution of the scaling limit. As for the ferromagnetic case, there are two particles 
A and A, but for the antiferromagnet the fundamental charged field which creates A is S 
(rather than a), and its odd sublattice parity forbids the fusion AA —^ A characteristic of the 
ferromagnetic solution (I76p . 

The above symmetry identifications do not depend on b, and this suggests that different 
values of b may describe the three-state Potts antiferromagnet on other lattices. It is known 
that at T = 0 the model is critical also on the Kagome lattice and allows for a Gaussian 
continuum limit |59] . 

5.5 Correlation functions 

The correlation functions are determined in the form ()16p if the form factors (|15l) are known. 
Remarkably, the form factors can also be determined when the S-matrix is known exactly. 
Indeed, in a disordered phase with factorized scattering the form factors of a scalar field local 
with respect to the order field satisfy the equations IMIIII] 

... ,0*,0*+i, ...,9n) = S{9i - 0,+i)F„^(0i,.. .,9i+i,9i ,... ,0„), (90) 

F^{9i + 2i7r, 92,..., 9n) = F^{92, ...,9n, 9i), (91) 

n 

ResB:=e+i-.F^+2{0’, 9,9i,...,9n) = i[l-W S{9 - 9k)]F^{9i, ...,9^), (92) 

k=l 

where we are considering considering the case of a theory with a single particle species (and 
then a single two-particle amplitude S{9)) which is sufficient to illustrate the general logic; we 
denote by F^ the matrix elements ()15p as functions of rapidities. Eq. (I90p simply states that in 

^^More generally, the critical temperature of the antiferromagnetic q-state Potts model on the square lattice 
was determined in m- 
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an integrable theory the exchange of the positions of two adjacent particles on the line produces 
a scattering amplitude. Eq. m is related to crossing and its meaning can be understood 
considering the case n = 2 . Then the equation can be written as ( 6 * 1 +f 7 r, @ 2 ) = F^( 02 ) 
and, considering that form factors of scalar fields depend on rapidity differences, says that Ff* 
takes the same value at i-K + 9 and zvr — 9, with 9 = 9i — 02 - Recalling Fig. [7] this means that, 
contrary to the scattering amplitude, the form factor does not exhibit the crossing cut in the 
s-plane, but only the unitarity cut associated to Eq. ([90]). This is understandable, since the 
scattering amplitude has four legs and the form factor only two. 

Eq. (|92l) prescribes the presence of a pole when charge, energy and momentum of two particles 
sum to zero. It is possible to set up a formal argument to see that these annihilation or 
kinematical poles are a peculiar two-dimensional implication of the general fact that matrix 
elements of with m particles on the left and n particles on the right reduce, through crossing, 
to plus “disconnected parts” due to annihilations [61] . We will show in section 17.11 that 

kinematical poles are needed on physical grounds, since they account for phase separation in 
two dimensions |62|. In theories exhibiting bound states form factors inherit from the scattering 
amplitudes the poles (|75l) and this allows one to write an additional residue equation relating 
F,^]^ to F^. Kinematical and bound state poles are the only singularities of the form factors as 
functions of rapidities. 

Here we illustrate the solution of the form factor equations for the Ising model {S{9) = — 1), 
which minimizes the technicalities remaining non-trivial. In particular, above Tc the order held 
cr{x) creates the elementary excitations, so that spin reversal symmetry requires that a couples 
only to an odd number of particles. All the descendants of a, however, share this property 
and the form factor equations (I90H92P necessarily allow for inhnitely many solutions for a hxed 
symmetry sector. It is a consequence of (fTTl) that, among the scalar descendants, the helds 
{dB)^(j) have form factors behaving as ^ ~ as \9i\ —>■ 00 . It is then natural to 

expect, and can actually be shown more generally [631164] . that the primary held 4> corresponds 
to the solution with mildest asymptotic behavior at high energies in the given symmetry sector. 
It is not difficult to check that for a such a solution is unique and reads |65] 

F2Vi(0i,---,02n+i) = *"Ffntanh^^, T > T,, (93) 

i<j 

where Ff is a constant by relativistic invariance. 

When we move to T < Tc we need to remember that, although S is still —1, the excitations 
are kinks interpolating between the degenerate vacua |0+) and |0_). The held a is topologically 
neutral and couples to states starting and ending on the same vacuum, i.e. with an even number 
of kinks. The fact that the kinks are non-local with respect to a (they are created by the disorder 
held /i of section [T3]) leads to a slight modihcation of Eqs. (1911) . (I92]) . Indeed, taking into account 
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vacuum indices (|MT) reads 

(O+|cT(O)|i^+_(01 + i7r)K.+ {e 2 )... i^-+(02i)) = 

{ 0 .\a{ 0 )\K.+ {e 2 )... K.+ {e2j)K+-{9i - ITT)) = 
-(O+|cT(O)|i^+_(02)...^+-(02,)^-+(0i -ivr)) (94) 

with the last equality following from spin reversal symmetry. Hence, if denotes form factors 
on the vacuum |0+), the r.h.s. of (I9ip acquires a minus sign. Moreover, when 6 i = 62 , the 
particle with rapidity 62 can directly annihilate that with rapidity 9i + iir in the hrst line of 
(1941) . while in the second line it has to scatter with 2j — 2 particles to reach and annihilate the 
particle with rapidity 9i — ivr, thus producing a factor These two different annihilation 

paths correspond to the two terms in the r.h.s. of (f9^ . so that the second one picks up the 
minus sign of the last line of (|94p . This is why the solution for a reads 

(01,..., 02n) = i"F^ n tanh , T < T, , (95) 

i<j 

where Fq = (0+|(t| 0+) is the spontaneous magnetization. The high energy limit, where the high 
and low temperature phases meet at the critical point, can also be used to show [66p [64] that 
T^(0i, 02 ) /^o 01 — 02 00 , a relation that fixes the relative normalization of form 

factors in the two phases. The energy density e, on the other hand, is invariant under spin 
reversal and topologically neutral below T^, it satisfies the same form factor equations above 
and below Tc with solution F^ oc 5 n ,2 sinh^^-^, consistently with the fact that it is a bilinear 
in the free fermions. 

The form factors (|93p and (I95p determine (cj(x)(t( 0)) in the Ising model through the large 
distance expansion (|16p . For this case the expansion allows a resummation and the result 
can be exactly expressed in terms of a solution of a differential equation of Painleve type, as 
originally shown in [67] taking the scaling limit of the lattice result. For theories with interacting 
particles, however, resummations are not known, even if all form factors are available. It is then 
very relevant for practical purposes to see how effective partial sums of (I16p can be. Taking the 
example of the susceptibility X = / d?x{a{x)cr{0)) in the Ising model, the leading contribution 
to the large distance expansion of the correlator comes from F[ above Tc and from below 
(the vacuum contribution is subtracted in the definition of the susceptibility). This two-particle 
approximation (i.e. taking into account contributions up to two particles) yields F+/F_ ~ 37.699 
for the universal ratio of the susceptibility amplitudes defined in (Hj), to be compared with the 
exact resuliP^ F_|_/F_ = 37.6936.. |67j . The origin of such a remarkable agreement has to be 
understood as follows. The truncated form factor sum for a correlator ($(3:)‘h(0)) becomes exact 
for |a;| large, but certainly fails to reproduce the exact divergence at short distance; 

the contribution of short distances, on the other hand, is suppressed upon integration over 
(Px = 27r|x|d|xl, and the suppression increases as A$ decreases. The order field, in particular, 

^^See |68l I50| for the complete list of exact amplitude ratios in the two-dimensional Ising universality class. 
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q 

field theory (2pa) 

lattice 

1 

160.2 [71] 

162.5 ±2 [72] 

2 

37.699 [69] 

37.6936.. [67] 

3 

13.85 [69] 

13.83(8) [7^174] 

4 

4.01 [69] 

3.9(1) m 


Table 2: Results for the universal susceptibility amplitude ratio r_|_/r_ in the g-state Potts 
model from field theory in the two-particle approximation (2pa) and from the lattice; q = 1 
corresponds to the amplitude ratio for the mean cluster size in percolation. 


has the smallest dimension and yields the best convergence for the integrated correlator. In this 
way the two-particle approximation becomes extremely effective in integrable theories, where 
two-particle form factors are not too difficult to compute, both for ordinary particles (see e.g. 
[50l[5l]; also m for a different approach) and for kinks |69[l7n| . In Table [2] we report the results 
of the two-particle approximation for r+/r_ in the g-state Potts model and compare them with 
the lattice results; for q ^ 2 the latter are obtained through series expansions or Monte Carlo 
simulations. 

In the g-state Potts model above Tc the probability that two sites have the same color a is 
given by 

('5s(xi),a'5s(x2),a) = -P2{xi,X 2) + - P2{xi,X2)) , (96) 

where P 2 is the probability that the two sites are in the same cluster. Recalling the definition 
([ 6 ]), we obtain 

{(Ja{xi)aa{x 2 )) = (g - 1)^2(a^i, X 2 ), (97) 

so that the susceptibility is proportional to the mean cluster size. When evaluating PU]) below 
Tc one has to take into account the possibility that the sites belong to the infinite cluster, but 
P7I) is recovered in the limit g —>• 1 relevant for percolation, with P 2 the probability that the 
sites are in the same finite cluster. Hence, with Ao-|q=i = 5/96 and A£|g=i = 5/8 from Table [H 
we obtain that the mean size of finite clusters in percolation diverges with the susceptibility 
exponent 7 = (1 — 2Ao-)/(l — A£)|q=i = 43/18 as p —)• Pc- (r+/r_)q=i, on the other hand, 
gives the universal amplitude ratio of the mean size of finite clusters on the two sides of the 
percolation transition; the agreement of the theoretical and numerical results for this number in 
Table [2] has to be regarded as particularly remarkable since its value had been controversial for 
decades, with numerical estimates spanning different orders of magnitude (see [771 EH]). 

6 Crossover behavior and massless scattering 

6.1 Dilute Ising model 

We saw in section 12.21 that for a given internal symmetry group G there may exist multicritical 
points characterized by the presence of more than one G-invariant relevant field, so that the 


31 







Figure 10: Schematic phase diagram of the Ising model with vacancies; the pure model corre¬ 
sponds to A = —oo. The first order (dashed) and second order portions of the phase transition 
line meet at a tricritical point. 


realization of such a multicriticality requires the tuning of more than one parameter. The basic 
example is obtained considering the Ising model and allowing for the presence of vacant sites, 
i.e. considering the Hamiltonian 

^diluteIsing — J ^ ^ O'iO'jtitj -\- A ^ ^ ti , CTj — il , ti — 0; 1 j (^8) 

ihj) * 

where ti = 0 introduce a vacancy at site i and A controls the vacancy density; as A ^ —oo 
the configurations with vacancies are suppressed and the undilute case is recovered. Dilution 
preserves spin reversal symmetry, which is observed to break spontaneously along a line in the 
(A)-plane, preserving the second order character of the undilute case up to a critical value 
of the vacancy density, and becoming first order above this value (Fig. llOjl . The point on the line 
where the order of the transition changes provides the simplest example of a tricritical point, 
i.e. a scale invariant point where a second symmetry-invariant field, the dilution e' conjugated 
to A — Ac, is relevant together with the energy density s conjugated to J — Jc- Adding e' to 
the action of the tricritical point induces a shift along the phase transition line, the first or the 
second order portion depending on the sign of the coupling. The correlation length is infinite 
along the second order portion of the line, but there is no scale invariance, since e' is relevant and 
its conjugated coupling dimensionful; moving down along the line, correlators cross-over from 
the power law behavior ruled by the scaling dimensions of the tricritical point to that ruled by 
the scaling dimensions of the ordinary Ising critical point. We call crossover line the transition 
line connecting two scale invariant points. 

In principle, the inhnite correlation length allows for a continuous, held theoretic description 
along the whole crossover line and its vicinity. As we discussed in section 13.31 in two dimensions 
the Ising tricritical point corresponds to the case t = 4 of the series of minimal conformal 
models; it has central charge c = 7/10, e = 4>i^2 and = 0i,3- Then the counting argument 
of section [Q says, in particular, that the transition line, corresponding to a shift in the 4 >i ^3 
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direction, is integrable in the continuum. The crossover portion, which has infinite correlation 
length, must correspond to an integrable theory with massless elementary excitations. 

Massless particles are normally excluded from treatments relying on the analytic properties 
of the ^-matrix m- The basic reason can be understood recalling Fig. [2l where as the mass 
m tends to zero infinitely many branch points collapse in the origin, thus making the analytic 
structure untreatable. In an integrable theory, however, the scattering amplitude has only 
two cuts in the s-plane (Fig. [7]) and the limit m —)• 0 can be taken in a controlled way [79]. 
Within the rapidity parameterization p = msinh/1 of the momenta, one can take m = Me““, 
a +00, and obtain right movers with e = p = (M/2)e® for (3 = 6 + a, and left movers with 
e = —p = (M/2)e“® for (5 = 9 — a, with M a mass scale. For the scattering of a right mover 
with rapidity 0 i and a left mover with rapidity 62 , the relativistic invariant (|lip becomes 

s = ^ (99) 

a relation which maps the branch point s = 0 at 9\2 ^ 61 — 62 = —00. This requires that the 
upper and lower edges of the unitarity cut, which in the massive case were mapped by (l69]l to 
real values of 612 with opposite sign, now correspond to the boundary values of two different 
functions, S{ 9 i 2 ) for the upper edge and 5(6*12) for the lower edge, with 6 real; the upper and 
lower edges of the crossing cut (s < 0) are obtained shifting 612 by ivr. With these notations 
and for a single particle species the massless limit of the unitarity and crossing equations (|13p . 
(IT^ reads 

,s(0)5(0) = i, s{e) = s{e + i7r), ( 100 ) 

which can be combined into 

S{9)Sie -i7r) = 1. (101) 

For the crossover line in the dilute Ising model, the mass scale M appearing in (|99p is 
associated to the coupling conjugated to the relevant field e', which is proportional to 
The high energy limit M^/s —)• 0 in which this coupling becomes negligible is that towards the 
tricritical point, while the critical point of the undilute model is approached as M^/s ^ 00. 
Since we know that the latter contains a single particle species, which is a free fermion, one 
looks for the simplest meromorphic solution of (|101l) giving —1 when 9 = —00. This solution 
reads 

^ (5 “ f) ’ 

and the thermodynamic Bethe ansatz can be used to confirm that it leads to the correct central 
charge 7/10 in the high energy limit |79] . The amplitude exhibits a pole at 9q = —f7r/2 (cor¬ 
responding to So = —iM^) in the unphysical strip lm.9 G (—Z7r,0). According to the general 
interpretation of second sheet poles this corresponds to a resonance, i.e. to an unstable 
particle with lifetime proportional to l/\/—Im sq = 1/M. As in the massive case, the scattering 
solution (11021) can be used to compute form factors and correlation functions |80j . 
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Figure 11: Phase diagram for two-dimensional Ising percolation. An infinite cluster of up spins 
exists above the percolation line. The percolation transition changes from first order (dashed) 
to second order at the Curie point. 


6.2 Correlated percolation 

In previous sections we considered the case of random percolation, i.e. the percolation problem 
in which each site (or bonc@) is occupied with probability p independently of the others. It 
is relevant, however, to consider also the case of correlated percolation, in which there is an 
interaction among the sites [5]. Once again the Ising model provides the basic example if we 
consider the spin clusters obtained connecting nearest neighboring sites with spin “up”, i.e. 
equal to +1. Indeed, let us consider the case in which a magnetic field H is present, i.e. we take 
the Hamiltonian ([3]) and add the term —H Uj. Then for H = +oo a cluster of up spins covers 
the whole lattice, while no up spins are left at H = —oo, indicating that a percolation transition 
takes place in between at some value Hc{l/J) (we introduce the notations J = J/T, H = H/T). 
We know that Hc{0) = 0, since at zero temperature the system is in the ferromagnetic ground 
state selected by the sign of H. For J = 0, on the other hand, sites are uncorrelated and the 
transition takes place when the probability that a spin is up coincides with the random site 
percolation threshold p^, i.e. Hc{oo) is determined by the equation /(2 cosh Hc{oo)) = p^. 

It is natural to expect that He is a monotonic function of the interaction among the spins, and 
it can be shown |81] that the existence of a spontaneous magnetization implies the presence of 
an infinite cluster. Since two-dimensional lattices have Pc > 1/2 [5], and then Hc{oo) > 0, these 
considerations lead to the phase diagram of Fig. [TTJ While for J > the percolation transition 
is, like the magnetic one, of the hrst order, the points = (j“^,0) and {oo, Hdoo)) are 

hxed points for percolation (correlated and random, respectively). The second order transition 
line connecting them is the crossover line; the linear size of clusters diverges along this line0, 
while the magnetic correlation length is everywhere finite away from the Curie point (jT^,0). 
This means that the field theory which describes percolation of spin clusters is not that which 
describes the local magnetization. 

In order to progress, let us recall what we saw for random percolation. In that case an 
auxiliary g-state Potts spin variable was associated to each lattice site, so that the “geometric” 

^^Site and bond percolation belong the same universality class. 

^^For the triangular lattice, which has p(l = 1/2, and then Hc{oo) = 0, the transition line remains at H = 0 for 
any J. 
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percolation transition was mapped onto the thermal Potts transition, in the limit g —>• 1 in which 
the auxiliary Potts degrees of freedom are removed. Exactly at g = 1 there are no degrees of 
freedom, but at g = 1 + e, e <C 1, we can extract the quantities relevant for percolation, e.g. the 
two-point connectivity which, recalling (1971) . is given by P 2 = A similar procedure 

can be used to track the Ising spin clusters, associating the auxiliary g-state Potts variable Sj 
to a site i in which the Ising spin is up. In practice one adds to the Ising Hamiltonian the 
term —JYl{ij) where = (uj -|- l)/2 = 0,1 [82l [83], and obtains in this way a dilute 

g-state Potts model, of which the dilute Ising model of the previous subsection is a particular 
case (g = 2). As for g = 2, the Potts model possesses, for real g < 4, a tricritical point and 
a crossover line connecting the tricritical point to the critical point of the undilute model [84] . 
For g = 1 -|- e this line describes the second order percolation transition line of Fig. [II][83|. 

Moreover, also for g 7 ^ 2 the held (^ 1^3 is responsible for the departure from the tricritical point 
towards the critical one cu, so that the whole transition surface spanned by g is integrable, in 
particular the percolation transition line at g = I -|- e [85] . The Sg-invariant massless scattering 
theory is not known for generic g, but the results of [86| are sufficient to show [85] that the 
amplitude Srl{0) possesses a resonance pole located at 9q = —i7r{t — 3)/2, where t determines 
the central charge c = 1 — 6/[t(t-|-1)] along the tricritical line; g = 2, with c = 7/10, corresponds 
to t = 4 and 6q = —i7rj2, as we saw in the previous subsection. On the other hand, at g = 1 
the only degrees of freedom left are the vacancies, which are the original Ising variables, so that 
we have c = 1/2 on the tricritical line, i.e. t = 3 and Oq = 0. The analytic mechanism behind 
the phase diagram of Fig. [ 11 ] is now clear. At g = 1 -|- e the crossover line corresponds to the 
presence of an unstable particle with lifetime proportional to 1 /e, plus a number proportional to 
e of massless particles (the Potts degrees of freedom); these massless particles yield an infinite 
connectivity length and account for the existence of the second order percolation line. Exactly 
at g = 1 we are left with the theory which describes the local magnetization properties of the 
Ising model; there are no massless particles left, while the former unstable particle has become 
stable and accounts for the hnite magnetic correlation length. 

7 Phase separation 

7.1 Order parameter profile 

So far we considered systems dehned on the whole plane. In this section we start the study 
of effects produced by the presence of boundaries, beginning with the case of an inhnitely long 
strip of with R. For discrete internal symmetry and T < Tc, a boundary field on the two edges 
of the strip can be used to drive the system in one of the degenerate phases as —>• 00 . We have 
already seen in section [5] that these phases are in correspondence with the degenerate vacua 
|0a), 0 , = 1,... ,n, of the associated quantum theory, and that the kinks Kab{9) provide the 
elementary excitations. On the other hand, if we place the strip as in Fig. [T2l and apply on 
both edges a boundary field favoring a phase a for x < 0 and one favoring a different phase b 
for X > 0, for R large we will have different pure phases at the far left and far right separated 
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Figure 12: Strip geometry. 


by a central interfacial region; R large means much larger than the correlation length ^ of the 
pure phases. Field theory provides a general and exact description of such a phase separation 
in the scaling limit [62l [87] . 

To see this we observe first of all that within the quantum picture with the vertical coordinate 
y as imaginary time, the boundary conditions play the role of initial and final states and can 
be expanded over the complete basis of asymptotic states of the bulk theory. In particular, the 
mixed boundary conditions specified above (we will refer to them as ab boundary conditions) 
correspond td^u 


\Bab{±R/2)) = 


I ^f{9)\Kab{e)) + ... 


(103) 


where we are first considering the basic case in which the vacua |0(j) and |0;,) are adjacent, i.e. 
there exists a single-kink excitation \Kab{0)) connecting them, and the dots stay for states with 
larger total mass yielding subleading corrections in the large R limit we are interested in. Then 
the partition function for the strip with ab boundary conditions can be written as 


Zab{R) = {Bab{R/2)\Bab{-R/2)) ~ ^|/(0)Pg—..iicosh0 _ 4== , (104) 

J 2-K ^I'KVTLabR 

where mab is the kink mass and we took the large R limit, as we will regularly do belov0. For 
uniform boundary conditions of type a the expansion of the boundary state \Ba) starts with 
the contribution of the vacuum |0a), and the corresponding partition function reads Za{R) = 
(B,(f)|B,(-|)) (0a|0a) = 1- Then we can determine the interfacial tension 

^ r 1 1 ^ab{R) 

T.ab = - hm - In = uiab ■ (105) 

R-^oo R Za\R) 


We now use the notations a{x, y) for the order field at the point (x, y), and {a)a = (0a|cr(a:, y)|0a) 
for the order parameter in the pure phase a. The order parameter profile in the middle of the 

^®In general, when the change from boundary condition a to boundary condition h takes place at xo -f 0, (ITO^ 
exhibits the factor 

^®Also for later use we recall the formula dxe~^^ = e'^ 
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strip is 


{a{x,0))ab = -^{BabiR/2)\a{x,0)\Bab{-R/2)) 

^ab 

~ j'^^{Kab{ei)\a{^MKab{e 2 ))e-^^^^+^^^^ (106) 

where O 12 ^61 — 62 and for simplicity we write m instead of niab- The matrix element inside 
the integral contains a disconnected term 2 Tr 6 { 6 i 2 ){o')b (corresponding to annihilation of the two 
particles) which contributes an additive constant to the profile; due to crossing, instead, the 
connected part is the two-particle form factor (06 |it( 0, O)|it'ba(0i + i'^)Kab{ 02 ))- For the Ising 
model, recalling (|95p . this is given by F 2 {0i + in, 82 ), and exhibits the annihilation pole at 
61 = 02 with residue 2i((T)+ = *[(< 7 )+ — (cr)-]. This pole and its residue owe nothing to the 
integrability of the Ising model. Indeed, integrability simplifies the scattering, but the matrix 
element in (jI06p is taken between states involving a single particle, and then no scattering. The 
pole then is always present as a consequence of the non-locality of the kinks with respect to the 
order parameter, and for a generic model the residue is i[{(y)a — (< 7 ) 6 ] [SS]- We write 


{Kab{0l)W{^-,^)\Kab{02)) connected 


,(0-)a - {(T)b 


0 


12 


+ E' 

n=0 


, an 
n "12 1 


(107) 


take the derivative of (|106l) with respect to x in order to get rid of the pole, perform the 
integrations over rapidities, integrate back over x with the boundary condition ((7(-|-oo, O))^;, = 
{(T)b, and obtaiiJ^ [62] 

{a{x, 0))ab = ^ [(C7)a + {a)b] - ^ [((7)a - {(T)b] +^0^^^g-2mxVR ^ ^ ^^gS) 

where the error function erf( 2 ;) = interpolates between —1 and 1, so that the 

profile interpolates between [a) a and {(j)b- Hence, the pole in (I107p is responsible for phase 
separation in two dimensions while, as we are going to see, the regular terms characterize the 

interface structure. For the Ising model (( 17 )+ = —{a)-, cq = 0) (jlOSp reduces to {a) _^ ~ 

(< 7 )+ eT:^{^J^^x), a result originally obtained from the lattice in [88] (see also [89] and references 
therein). 


7.2 Interfaces 

The notion of interface emerges from (11081) in the following way. In first approximation we 
think of it as a simple curve connecting the boundary condition changing points (0,±i?/2), 
intersecting the axis y = 0 only once and sharply separating two pure phases IFig. fTHhf : then we 
add corrections to this picture allowing for a structure localized on the curve. Hence, we write 
as 

aab{x\u) = 0(u - x){a)a + 0(x - u){a)b -|- Ao 6 {x - u) + Ai 6 \x - u) + ... (109) 

take into account that, for phases a and b playing a symmetric role, f{9) = /(O) -I- 0{6^). 
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(a) (b) (c) (d) 

Figure 13: Configurations of a single interface (a,b,c) and a double interface (d). 


the order parameter at the point (x, 0) in a configuration in which the curve intersects the central 
axis at (tt,0); 0(x) is the step function equal to 1 if x > 0 and to 0 if x < 0, while the delta 
functions account for the structure on the curve. Calling p{u)du the probability that the curve 
passes in the interval (rt, u + du) on the axis y = 0, we then have 

{a{x,0))ab = j duaabix\u)p{u), (110) 

which coincides with the field theoretical result (I1U8|) for p{x) = 

We see that the term proportional to cq in (jlOSp accounts for the first deviation from sharp 
phase separation, an effect which is natural to ascribe to bifurcation and recombination of the 
interface (Fig. I13bi: it requires three different phases, consistently with the fact that cq = 0 in 
the Ising model. For the latter, the leading correction will come from trifurcations (Fig. fT3b). In 
the three-state Potts model, on the other hand, {( 73)1 = (< 73)2 and cq / 0 [62], so that the leading 
variation in the profile ((73(x,0))i2 is produced by the Gaussian term in (11081) . corresponding to 
small drops of phase 3 forming on the interface between phases 1 and 2; the effect is suppressed 
as as i? —)■ 00 . 

Using (1141) the order parameter calculations that we performed for y = 0 can be repeated 
for any fixed \y\ < R/2, in the limit i? —)• 00 |87|. From the order parameter one then obtains 
the probability density 

, . 1 12m 

for the passage in the interval (x, x -|- dx) on a line of constant y; we introduced the definitions 

K^^l-4yyR\ (112) 

so that (|llip exhibits a y-dependent variance which shrinks to zero as y approaches the boundary 
condition changing points, as it should. 

The above results have been obtained for adjacent vacua |0a) and |0b). In some cases, 
however, connecting the two vacua requires a n-kink excitation with n > I. The g-state Potts 
model can be used to illustrate this case. For g < 4 the transition is of the second order and 
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Figure 14: Vacuum and kink structure for the three-state Potts model below (a), the dilute 
three-state Potts model at the first order transition (b), and the Ashkin-Teller model below Tc] 
in the latter case the dashed lines correspond to kinks present only in the attractive regime. 


below Tc there are q degenerate vacua which, by symmetry, are adjacent to each other (Fig. [TTki. 
In presence of vacancies for (7 < 4 there is also a tricritical point and a first order transition line 
(Fig. [T 0 ]l on which the ferromagnetic vacua | 0 j), i = 1,... ,q, are degenerate with the disordered 
vacuum |Oo). In this case the kinks KiQ and Kqi are the elementary ones, since two ferromagnetic 
vacua can be connected by a two-kink excitation KioKQj. This first order transition is integrable 
in the scaling limit and the exact S-matrix (DHI shows that KioK^j do not form bound states, 
so that the ferromagnetic vacua are not adjacent to each other (Fig. 114b). The same conclusion 
holds for the pure model with T = Tc and g = 5, which also is a first order transition point. 
Strictly speaking, the scaling limit for g > 4 is only defined in the limit g —)■ 4+ in which the 
correlation length tends to infinity, but for g = 5 (and actually up to g ~ 10) the correlation 
length is much larger than lattice spacing and an effective S'-matrix description is possible and 
accurate m- 

When the lightest excitation connecting |0a) and |0fe) is a two-kink one, the boundary states 
(|103[) become 


Bcb{±R/2)) = ^ 


^ J deide2facbi0l,e2) \Kac{ei)Kcb{92)) +... 


(113) 


and give rise to the double interface of Fig. [T3l i. The calculation of the order parameter profile 
starting from these boundary states is a generalization of that performed above for adjacent 
phases [57]. For the Potts first order transition points discussed a moment ago one obtains, in 
particular. 


(cj3(x,0))i2 ~(^1 - erf(z)e , z=\J^x (114) 

as the leading term for mR S> 1. Now limij^oo(c’'3(0,0))i2 / (< 73 ) 1 , in contrast to what we saw 
for the pure Potts below T^; the difference is that now a macroscopic drop of the disordered 
phase forms between the two interfaces. 

As for the single interface, the results for the order parameter profile can be extended to 
y / 0 and used to obtain the probability density that the two interfaces pass at points xi and 
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X 2 on a line of constant y [871. This reads 


p{xi,X2,y) = e (115) 

the quadratic prefactor prevents the interfaces from passing at the same point and provides the 
effective repulsion which inflates the bubble of the third phase. In the field theory derivation it 
arises quite generally from the fermionic low-energy behavior of the scattering amplitudes m, 
and this explains the appearance of (IllSp in the framework of the lattice Ising model, suitably 
reinterpreted to account for the lack of a third phase |9Tl [92] ; (11151) is also at the basis of the 
heuristic treatments of mutually avoiding interfaces [93] . 

Let us remark, before closing this subsection, that the interfaces we are concerned with 
here and below are associated to separation of coexisting ordered phases and that, as we have 
shown, their properties can be derived directly in the continuum from the determination of the 
order parameter. The situation is different if the system is exactly at a second order phase 
transition point, where there are no coexisting phases and (cr(x, O))^;, goes to zero for any x as 
i? —>■ oo. Still, if we consider the Ising model at on the triangular lattice with boundary 
conditions, any configuration exhibits a path (on the dual honeycomb lattice) with positive 
spins on its left and negative spins on its right connecting the boundary condition changing 
points (0,±i2/2). It has been shown in recent years that the passage probability for such a 
path can be determined exactly within the mathematical theory of stochastic curves subject to 
the constraint of conformal invariance [94] . In conformal held theory the same result emerges 
[95] as the solution of the second order differential equation associated to the degenerate held 
which induces the change of boundary conditions at the points ( 0 , ±i?/ 2 ) (see e.g. [96l [97] 
for reviews). When passing from Ising to the three-state Potts model, the paths connecting 
the boundary condition changing points can branch even on the honeycomb lattice, and this 
case has been studied at criticality using Monte Carlo [95] and transfer matrix [99] techniques. 
Looking for a relation between these lattice paths at criticality and the off-critical interfaces 
characterized earlier in this section directly in the continuum appears as an interesting subject 
for future research. 

7.3 Wetting 

We have seen how different vacuum connectivities give rise to two different regimes of phase 
separation: a single interface with small bubbles along it (adjacent vacua) or a double interface 
enclosing a macroscopic layer of a third phase (non-adjacent vacua). The change of a parameter 
of the system may induce a transition from the first to the second regime which, thinking of 
liquid droplets which condensate to form a layer, is known as wetting transition (see e.g. m). 

The Ashkin-Teller model provides an exact illustration of such a wetting transition [57] . The 
lattice Hamiltonian (|89p shows that the model corresponds to two Ising models coupled through 
the product of their energy densities, and we have seen that this gives rise to a line of critical 
points in the space of the couplings J and J 4 , and that the scaling region around this line is 
described by the sine-Gordon model. In the regime in which the reversal symmetries of both 
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Ising spins are spontaneously broken the model possesses four degenerate ferromagnetic vacua 
| 0 (ji,a 2 )! fundamental kinks interpolating between them are associated to the soliton 

A = Ai + iA2 and antisoliton A = Ai— iA2 of sine-Gordon, since Aj corresponds to the kink of 
the j-th Ising model. Sine-Gordon possesses the parameter (I87p which serves as the coordinate 
along the critical line, with the free fermion point ^ = vr corresponding to the Ashkin-Teller 
decoupling point J 4 = 0. In the sine-Gordon repulsive regime ^ > vr (J 4 < 0) Ai and A 2 are 
the only kinks and connect the vacua along the sides of the square in Fig. [T4h : vacua separated 
by the diagonal are not-adjacent since connecting them requires a two-kink state A 1 A 2 . In 
the attractive regime S, < tt (J 4 > 0), instead, A 1 A 2 form bound states and all vacua become 

adjacent. In particular, if we consider the strip with boundary conditions-for x < 0 and 

-|—I- for X > 0, we have a single interface with interfacial tension S_= mi = 27T7-sin(^/2) 

(the lightest bound state mass given by ( 1881) 1 in the attractive regime, and a double interface 
(enclosing a mixture of the phases and —h) with interfacial tension 2 m in the repulsive 
regime. There is then a wetting transition at J 4 = 0 corresponding the unbinding of the bound 
state (mi = 2 m). 

A different mechanism, kink confinement [1011 1102] . is responsible for the change of the 
interfacial structure when we start from a first order transition point in which ordered phases 
coexist with the disordered one (we gave above the examples of the dilute {q < 4) and pure 
{q > 4) Potts model), and then move the temperature below Tc. Let us first discuss kink 
confinement independently from wetting, within the simplest context of the Ising model. If 
we start at T < Tc and add to the action a symmetry breaking term hj d?xa{x), with h 
small, this removes the degeneracy of the two vacua inducing a difference S£ ~ 2/i(cj)+ in their 

energy densities. An immediate consequence is that the kink _is no longer a stable particle. 

Technically this appears through the fact that for small h the correction to the kink mass is 

given by h{K^ _(0)|(T(O)|iL_= hF^{6 + iir, 6), which is infinite due to the annihilation pole 

in (|95p : as expected, explicit symmetry breaking removes topologically charged states from the 
spectrum and, as a consequence, forbids phase separation at equilibrium. On the other hand, if 
the sign of h is such to favor | 0 +) as the new ground state, the fate of the topologically neutral 

composite iL_|_ K _^ is different. Indeed, to first order in h the energy of a configuration in which 

the two particles are separated by a distance R on the line grows with R as R5£. This is a 
confining potential which prevents the kinks from moving infinitely far apart and gives rise to a 
string of “mesons” (kink-antikink composites). The lightest of these correspond to confinement 
of slow kinks and their masses are determined by the Schrodinger equation with the linear 
potential as 

Mn ~ 2 m H-— , (116) 

where n = 1 , 2 ,..., are positive numbers determined by the zeros of the Airy function, 
A\{—Zn) = 0. The spectrum (jll 6 p was first obtained from the study of the spin-spin correlation 
function on the lattice [103] : the evolution of the mesonic spectrum as the magnetic field becomes 
stronger can also be followed [lUll 110411105| . In the three-state Potts model a symmetry breaking 
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field favoring the vacuum |0i) leads also to the formation of “baryons” from the conhnement of 
three-kink states Ki 2 K 2 sK^i [TUmfWlfTUg] . 

Going back to the q-state Potts model at a point of coexistence of the ferromagnetic phases 
i = 1,... ,q with the disordered phase 0, we know from Fig. 114b that boundary conditions on 
the strip of type i for x < 0 and j for x > 0 lead to the formation of a wetting intermediate 
disordered phase produced by the state Ki^Koj. When moving the temperature slightly below 
Tc, the energy density of the disordered vacuum exceeds that of the ferromagnetic vacua by 
5E oc (Tc — T), and the state KioK^j gives rise to stable “mesonic kinks” with the masses 
(11161) . Thus we go back to the connectivity depicted in Fig. [T4b . i.e. the phases i and j are 
separated by a single interface. 

7.4 Quenched disorder 

The results seen so far allow few considerations about quenched disorder which, in a slightly 
different form, were originally developed in |109j . Ferromagnetic quenched disorder corresponds 
to replacing the constant coupling J among sites in the lattice Hamiltonian with site-dependent 
couplings Jij randomly distributed (e.g. in a Gaussian way) around a positive mean; one then 
write the partition function Z{{Jij}) summing over spin configurations for a given realization 
of the disorder, and takes as free energy of the system the average over the disorder of — In Z. 
A usual way to deal with such an average is to writ 

_ Yn _ 1 

\nZ = lim-, (117) 

n^O n 

so that one is led to consider the partition function Z” of n copies (replicas) of the system with 
the same bond configuration. For weak disorder one can then argue that the average results into 
an effective interaction among n replicas of the pure system {Jij = J) (see e.g. [7]). Apart from 
the averaged free energy, the quantities Z^ = determine the fluctuations of the random 

variable In Z. 

Gonsider now an Ising model on the strip geometry of Fig. 1121 with boundary conditions — for 
X < 0 and -|- for x > 0, in presence of weak quenched disorder. Following the replica treatment, 
for temperatures below the critical temperature of the pure system and for R much larger than 
the correlation length, we find n interfaces connecting the boundary condition changing points 
(0,±ii/2). Assuming that, as indicated by perturbative calculations [llOj . the phase transition 
persists in the disordered system, we have to think that the interaction among the n interfaces 
induced by the disorder is attractive and binds them into the single interface required for a 
system with two degenerate vacua; for weak disorder the binding will be weak. It is well known 
(see e.g. m and references therein), that in 1 -|-1 dimensions the binding of n identical particles 
becomes model-independent for weak attraction, since it approaches the non-relativistic limit 
of pairwise delta function interaction. Hence, in particular, we can read off the result from 
the sine-Gordon model of section 15.31 where we saw that the particles Bn, with masses 

^®We denote by X the average of X over the disorder. 
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(a) 


(b) 


Figure 15: Reflection in the direct channel (a) is related to pair emission in the crossed channel 
(b). 


given by (j88p . are bound states of n particles Bi. In the weak binding limit ^ 0 we obtain 

uin ^ M{n — n^^^/24) for the interfacial tension in the disordered system. Recalling (jl04p we 
then obtain Z'^_^ oc for R large. The fact that the leading effect of the disorder 

comes with the term in the exponent is a relevant result for two-dimensional systems with 

quenched disorder, since it can be shown that it leads to characteristic fluctuations of In Z _of 

the form (rather than at large R [109] . 


8 Interfaces at boundaries 

8.1 Boundary scattering 

Critical and near-critical phenomena at boundaries provide an important chapter of the theory of 
universal behavior in statistical physics (see e.g. [11211113] ). In two dimensions, the critical case 
is described by boundary conformal field theory [114] . which enables one to obtain exact results 
for boundaries of any shape exploiting conformal mappings and infinite dimensional conformal 
symmetry. Here we are interested in the near-critical case, and in this subsection consider the 
interaction of the bulk degrees of freedom with the boundary, for the basic case of the half¬ 
plane with coordinates {x > 0, y) and translation invariant boundary conditions. Within the 
scattering framework in which y is the imaginary time, this corresponds to consider particles 
moving on the half-line a; > 0; they scatter among them but also with the boundary, which one 
can conveniently think as an inhnitely massive particle stuck at x = 0. It is possible to follow 
a logic similar to that used for the scattering in the whole two-dimensional space-time and to 
look for exact solutions characterized by complete elasticity and factorization of the scattering 
amplitudes [115] fsee also m)- As for the whole plane, these solutions exist and describe the 
main interesting models; they also explicitly exhibit features, such as boundary bound states, 
which are general and will be recalled at some points later. 

Integrability on the half-plane requires that on the whole plane to start with, and two-particle 
scattering is specified by a known amplitude S{ 6 i — 62 ) of the bulk theory; we are referring to 
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the case of a single particle species which is sufficient to illustrate the main ideas. Due to 
complete elasticity and translation invariance in the time direction, a particle traveling towards 
the boundary with rapidity 0 < 0 will bounce back with rapidity —9 (Fig. [TShl. and we denote 
by R{6) the amplitude for such a reflection process. Recalling that the boundary can be seen as 
a particle with zero rapidity, the unitarity equation 


R{e)R{-9) = 1 (118) 

for the reflection amplitude is a particular case of (I70p . The boundary version of crossing 
symmetry, on the other hand, is more subtle and can be understood as follows |115| . If we 
look at the process of Fig. [15^ from the “crossed channel”, keeping imaginary time along the 
y direction, it will appear as in Fig. 115b . i.e. as the emission from the boundary of a pair of 
particles with rapidities 0 and —9 (momentum, rather than energy, is conserved in the crossed 
channel). Hence, the reflection amplitude in the direct channel is related by analytic continuation 
to the pair emission amplitude, that we denote K{9), in the crossed channel. Since in passing 
from Fig. [T5h to Fig. 115b energy and momentum have interchanged their role, the analytic 
continuation following from (j68p is 

IIT 

K{9) = R{--9). (119) 

On the other hand, interchanging the position of the particles with rapidities 9 and —9 in the 
crossed channel produces a factor S{29)', then we have the relation 


K{9) = S{29)K{-9). 


( 120 ) 


which through (jllOp provides a second equation for R{9) |115j . 

For the Ising model, with S{9) = —1, the equations 
solution 


( 27r 

T ~ 2 


9\ K — i sinh 9 


admit the one-parameter 


( 121 ) 


K + i sinh 9 

corresponding, for T < Tc, to the presence of a boundary magnetic field h related to k as 


K = 1 -. 

2 m 


( 122 ) 


In this case the identification of the boundary condition is allowed by the fact that the boundary 
operator conjugated to h has dimension A = 1/2 |117j . and then is linearly expressible in terms 
of the free fermions of the bulk theory; this allows one to obtain (11211) and (11221) through free field 
techniques |115j . Notice that Ro{9) = —zcoth(^ —|) possesses a pole at 9 = i7r/2 which is absent 
for /i 7 ^ 0 and turns out to be a general feature distinguishing symmetry-preserving (or free) 
from symmetry-breaking boundary conditions. Together with minimality this normally allows 
one to identify solutions of the boundary scattering problem corresponding to fixed {h = oo) and 
free boundary conditions below (see e.g. |118j for the Potts model). Identifications below Tc 
are made easier by the fact that the bulk excitations are kinks whose trajectory divides the half¬ 
plane in Fig. [15] into magnetized domains which have to be compatible with the magnetization 
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state of the boundary. Information about the system above can usually be obtained by duality 
considerations (see m and below). 

Another instructive feature of (I12ip is the presence of a pole at 6 = iu, with u satisfying 
sinu = 1 ~ For < 2 m this pole is in the physical strip and corresponds to a boundary 
bound state, i.e. to a bound state formed by the boundary and the bulk particle which has to 
be interpreted as a stable excited state of the boundary; the difference between the energy of 
the excited boundary and that of the ordinary boundary is given by the energy of the particle 
evaluated at 9 = iu, i.e. by 

eB'— gb = mcosu . (123) 


8.2 Interface in a wedge 

Interfacial phenomena at boundaries are extensively studied for their relevance in applications 
involving adsorption of fluids on substrates (see e.g. [1201112111122j l. In this and in the next 
subsection we show how field theory provides general and exact results for universal properties 
in the two-dimensional case [12311124j . We start from the half-plane with the boundary located 
at X = 0 and consider the scaling limit below Tc in which two or more phases coexist; as usual 
the kinks Kab provide the elementary excitations in the particle description. We denote by Ba 
a boundary condition corresponding to a boundary field favoring phase a in the bulk. We can 
switch from boundary condition Ba for y > uq to boundary condition Bb for y < yo inserting at 
{x,y) = (0,yo) a boundary condition changing field yiab with the following dehning properties. 
First of all it has to couple to kink states connecting the vacua |Oa)o and |06)o, where the 
subscript 0 indicates that we consider the quantum system dehned for x > 0. The simplest 
matrix element of yab is 

o{>^aM^,y)\Kba{e))o = , (124) 

where J-q (6 ) gives the amplitude for the emission of a kink from the boundary condition changing 
point. For the time being we want to exclude the possibility that the kink remains stuck on the 
boundary, and then we require T'q (0) = 0 or, on general analyticity grounds, Bq{9) = c9+0{9^), 
with c a constant. 

For reasons to become immediately clear we observe, recalling (1681) . that a rotation which 
brings the boundary to form an angle a with the vertical axis corresponds in real time to a 
relativistic transformation which shifts rapidities by —ia, so that the kink emission amplitude 
Ta in the rotated frame is related to that in the original frame as 

= (125) 

in particular 

F^{e) r. c{e + ia), |0 |,|a|<l. (126) 

With this information we can consider, instead of the half-plane, the more general wedge ge¬ 
ometry of Fig. [161 The points (0, ±i?/2) are boundary condition changing points, such that a 
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Figure 16: Wedge geometry with boundary condition changing points at (0,±i2/2) and the 
interface running between them. 


boundary field favors phase a (resp. h) for \y\ > R/2 (resp. \y\ < R/2) on the wedge. With 
these boundary conditions the order parameter in the wedge, that we denote by {a{x,y))wahai 
tends for X —>■ oo to {a)a if R is hnite, and to {a)i, if R is inhnite. For mR large we then expect 
an interface running between the points (0,±ii/2), separating an inner phase b from an outer 
phase a, and whose average mid-point distance from the wedge diverges with R. 

The theory confirms these expectations as follows. The order parameter for \y\ < R/2 reads 


{cr{x,y))waba 


a(0a|/iab(0, f )cr(x, 

^Waba 



(127) 


where the subscripts ±q: indicate the different rotations performed for positive and negative y, 
and 


^Waba = a(0a|lta6(0, 7?/2)/X6a(0,-i?/2)|0a)-c 


f 


de 






g2 g mR 


(1 + mRa^ 


(128) 


2^- av-.- -avw - 2V^(mii)3/2 

here we considered the limit mR large which projects on the lightest (single-kink) intermediate 
state, and a small in order to use (I126p . Similarly we have 




^—mR r+oo 


L 




Zw^ba J-cx> (2vr) 

X Rf/M {Kab{0i)\a{0,0)\Ki,a{e2)) > 


(129) 


where we evaluate the order parameter field on bulk states implying that the boundary condition 
changing fields account for the leading boundary effects at large R. Recalling (I107p and (I112p 
we obtain 


dx{(^{x,y))Wa,, _ . (rn\I + {ayf 2 

{o')a-{cr)b + mRa"^) 

and, integrating back over x with the condition {a{+oo,y))waba ~ {^)a-, 


(130) 


{(T{x,y))w,t, ~ {a)b -b [(o-)a 




erf(x) 


x +V2^f 

^/7^ 1 -b mRa^ 


(131) 
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Figure 17: Passage probability density p{x\ y)/m for the interface from (jl30jl with mR = 25, a = 
0.04. The leftmost contour line corresponds to p{x] y) = 0, and then to the wedge. As expected, 
the passage probability density diverges as the pinning points (0, ±i?/2) are approached. 

for a = y = 0 and (a)a = —{(^)b this result coincides with that obtained in |125j from the exact 
lattice solution of the Ising model on the half plane. Proceeding as for phase separation in the 
strip, we can write 

rx roo 

{a{x,y))w,,ha ^ {^)a dvp{v;y) + {a)b dvp{v,y), (132) 

J X J X 

where p{v, y) is the probability that the interface passes in the interval [v, v + dv) on the line 
of constant ordinate y, which intersects the wedge at x = x{y)] it follows that p{x;y) coincides 
with (11301) . The requirement dxp(x;y) « 1 leads to VmRa <C 1. Notice that (I13np shows 
that p{x', y) vanishes for |y| = § + "f, which for the present case of small a are the coordinates 
of the wedge {x > —i?a/2); hence, the properties (|125l) . (|126l) that we identified in momentum 
space indeed lead to an impenetrable wedge in coordinate space. A plot of p{x\y) is shown in 
Fig. [m 

8.3 Boundary wetting and wedge filling 

In deriving the results (jl30p . (|131l) we considered a regime in which the kink (i.e. the interface) 
does not stick to the boundary. On the other hand, considering first the case of the half-plane 
(a = 0), we saw in section ISTTl that in some range of the boundary parameters the particle may 
form with the boundary a stable bound state with energy e^/ specified by (I123p . i.e. smaller 
than the energy eb + m cosh 9 of the unbound boundary-particle system. In the present case 
the particle is a kink and such a boundary bound state, let us denote it |0(j)o, corresponds 
to the phase b forming a thin layer adsorbed on the boundary (Fig. I18p . When this state is 
present it dominates the spectral decomposition of (11271) at large mR and produces an order 
parameter equal to (a) a for mx 1, no matter how large R is m- We know from (|105p 
that m is the interfacial tension between the phases a and b; on the other hand, es corresponds 
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Figure 18: When the kink rapidity approaches the resonance value iu the kink and the boundary 
form a boundary bound state in which phase b is adsorbed on the boundary as a thin layer; u 
becomes the contact angle of the phenomenological wetting theory. 

to the interfacial tension between the boundary and phase 6, and that between the excited 
boundary and phase a (Fig. I18|) . It follows that (|123p amounts to the condition which in the 
phenomenological theory of boundary wetting (see e.g. [122j ) determines the equilibrium of a 
drop of phase b between the boundary and phase a. Moreover, the bound state parameter u is 
identified with the contact angle of the phenomenological theory, i.e. with the angle that the drop 
forms with the boundary at the contact point. The wetting transition is phenomenologically 
identified with the vanishing of the contact angle, when the drop spreads on the whole boundary; 
in field theory it corresponds to the unbinding of the kink from the boundary which takes place 
at u = 0 |123j . 

For fixed values of the boundary parameters, u is a function of the temperature, which is 
related to the bulk correlation length ^ oc m~^ oc (Tc — T)“^, so that the condition u{To) = 0 
defines the wetting transition temperature Tq < Tc for the flat boundary. For the Ising model 
with a boundary magnetic field h discussed in section [8T] we have 1 — sinu = ^ = 
the last equality follows from vising = 1 and holds in the scaling limit. 

In order to consider the case of the wedge (a > 0), we recall first of all that matrix elements 
of helds on particle states inherit from the S'-matrix the poles corresponding to bound states. 
For the matrix element (|124p the pole induced by the boundary bound state takes the form 

iV^' 

= Q-in. (133) 

where is the kink-boundary-excited-boundary coupling. It then follows from (11251) that 
the boundary bound state pole of Fa{9) is located at 0 = i{u — a). Since the kink energy 
is always m cosh. 9 and 9 = 0 remains the scattering threshold, the wedge wetting (or filling) 
transition temperature r„ is determined by the condition 

u{Ta) = a . (134) 

For u < a the pole is located at Im0 < 0, namely outside the physical strip Im0 £ (0, vr) allowed 
for bound states; in such a case the kink is unbounded and phase b fills the wedge. The condition 
(jl34p . which follows here from the exact field theoretical analysis, is known phenomenologically 
|126j and provides the basic illustration of how the adsorption properties are affected by the 
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Figure 19: A cluster configuration on the rectangle including a vertically crossing cluster. 


geometry of the substrate |127l I128j . It relates the filling transition temperature in the wedge to 
the contact angle u{T) for the fiat boundary, a connection known as wedge covariance |129U130] : 
field theory shows that wedge covariance originates from the relativistic covariance expressed by 

(M HM]. 

9 Off-critical crossing clusters 

Let us consider random percolation on a rectangle of sides R and L (Fig. [19]). We call vertically 
crossing cluster a cluster touching both horizontal sides of the rectangle, and vertical crossing 
probability the probability that at least one such cluster is present. In the scaling limit 
is a universal function of R/^ and L/^, where ^ is the bulk correlation length, proportional to 
the average linear cluster size. Crossing probabilities have been originally studied at the critical 
point (occupation probability p = pc) |131] . where ^ = oo and Py becomes a function of a single 
variable, the aspect ratio R/L. For this case Cardy was able to determine Py exactly mapping 
it to a four-point function of boundary condition changing fields associated to the corners of the 
rectangle; since these fields are degenerate, in the sense explained in section 13.11 Py becomes 
the solution of a differential equation of hypergeometric type |132j , a result which later has also 
been obtained through mathematically rigorous methods |133j . 

Here we discuss the near-critical case and show how, building on ideas introduced in the 
last sections, Py can be determined exactly in the limit in which the sides of the rectangle 
are much larger than ^ m- As we regularly did to deal with percolation, we start with a 
mapping to the g —)• 1 Potts model. We recall that the factor q to the number of clusters in 
the cluster decomposition ([7|) of the Potts partition function is due to the fact that each cluster 
comes in q colors. When putting the model on a rectangle, we have to take into account the 
effect of boundary conditions. Let us call the Potts partition function on the rectangle with 
boundary conditions u, d, I and r on the upper, lower, left and right boundary, respectively; we 
will then denote by / the free boundary condition, and by a the boundary condition in which 
the spins are fixed to a color a. Then, for a ^ b, does not receive a contribution from 
configurations containing vertically crossing clusters. On the other hand, as <7 —>• I the weight 
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of a configuration becomes its probability, and one has 

= ( 135 ) 

a relation already used in m for the study of the critical case. In order to express ^ud away 
from criticality we consider first the limit L —>■ oo and use the boundary state technique of 
section rm We then write 


Z lr 
ud 


= {B^ 


-RH 


\Bd) 


Lr 


L 


oo , 


(136) 


where now {Bd^u) are boundary states corresponding to uniform boundary conditions (fixed or 
free), and the vertical boundary conditions at infinity have the role of selecting the states which 
can propagate between the horizontal boundaries. Consider the case T < in which there are 
q degenerate vacua |0a) and the bulk elementary excitations are the kinks Kab with mass m. 
Then fixed and free boundary states expand aa^u 


\Ba) = 


|0a) + I ^f{e)\Kab{-e)Kbam 


+ 


m = E 


|0a) + S E l^ab(O)) + - - - 

b^a 


(137) 

(138) 


respectively, where, due to translation invariance, all the states appearing in the expansion have 
zero total momentum, as originally illustrated in |115j . Then we have, in particular. 


[Zg]+ = {Bf\e-^^\Bf)X, = fmLe-^^ + 0(e-2™«), (139) 


where the vertical boundary conditions a ^ h forbid vacuum contributions, and we used {Kab{^)\Kbc{^)) 
2Tr6acS{0) = SacmL; we see how the L-dependence is reintroduced through the volume regular¬ 
ization of the delta function. The superscript -|- in (|139p recalls that the quantity is computed 
in the scaling limit below in the Potts model, which corresponds to the scaling limit above 
Pc in percolation; we will use the superscript — for the scaling limit below pc- It is not difficult 
to see that there is a one-to-one correspondence between lattice Potts configurations with verti¬ 
cal crossings and configurations without horizontal crossings on the dual lattice [3]; then Potts 
high-low temperature duality leads to m 


P- = 1-P+ = [Zff]+^,^AmLe-^^, (140) 

P+ = 1-P^ = l-P-\R^L-^l-AmRe-^^. (141) 


where Ph is the horizontal crossing probability, A = g‘^\q=i, and the results hold for mR and 
niL both much larger than 1. Recalling (I97p we have m = Ij^ below pc and m = 1/2^ above 
Pc, with ^ oc. \p — Pc\~^^^ the correlation length dehned by the exponential decay of the bulk 
two-point connectivity. 

^®See |134| for the expansion of boundary states in the basis of high-temperature excitations. 
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Recall now Fig. 115b . where a pair of particles with rapidities 6 and —9 is emitted from the 
boundary with emission amplitude K{0). If the rapidity difference 29 coincides with a value iu 
for which the two particles can form a bound state, K will have a pole of the form ^ [ns], 
where F is the three-particle coupling (Fig. [8]) and g the boundary-particle coupling. Clearly the 
latter is related to g in (11381) |115j . the precise relation being g = gl2 |135l I136j . On the other 
hand, the Potts model on the half-plane with fixed or free boundary conditions is integrable in 
the scattering framework of section 18.11 |118j . and this leads to nsi 

A = . (142) 

Notice, recalling ([7|), that lim^^i In zJJ gives the mean number of clusters. The same 
operation on z// gives the mean number of clusters which do not touch the horizontal sides of 
the rectangle, since the others are forced to have color a and to not contribute a factor q to the 
partition function |137j . Similarly, one can count also the clusters touching only the upper or 
lower side, and arrives to the expression 

yffyff 

- ^ ^ f f ^CLCL 

= lim dp In —jj (143) 

af fa 

for the mean number of vertically crossing clusters. When using (jl37p and (jl38p to evaluate 
the partition functions above pc in the limit of large L and R, the finite size correctionj^ are 
dominated by the single-kink contribution coming from zjj, and we obtain [119j 

~ lim dp In zll ~ lim dg ln[g -|- q{q — l)mLg^e~^^] ~ 1 -|- A mL (144) 

the first term is produced by the vacuum contributions and shows that above pc there is a single 
cluster which becomes infinite in the limit i7 —)• oo. Below pc, on the other hand, vertically cross¬ 
ing clusters become rare as R becomes large, and we have P~ = Froh{Np > 0) ~ Froh{Np = 
l)~iVT. 

Next-to-leading corrections are computed in m. where the comparison with numerical 
data can also be found. In principle, integrability allows an exact evaluation of partition func¬ 
tions on strips also for R not large, through the boundary version |138| of the thermodynamic 
Bethe ansatz; at present, however, the thermodynamic Bethe ansatz for the Potts model with 
continuous values of q is unavailable even in absence of boundaries (see |139j ). 

10 Conclusion 

In this article we illustrated how field theory leads to the exact description of universal prop¬ 
erties of two-dimensional statistical systems. We tried to convey in simple terms the origin of 
the specific effectiveness of held theory in two dimensions, and showed its wide applicability, 

^°Notice that (I136II determines the partition function up to a bulk contribution which does not affect our results 
for crossing clusters. 
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which allows one to treat within the same theoretical framework problems such as percolation 
and wetting, usually considered as uncorrelated. For the critical case, we stressed the role of 
conformal field theory, but also pointed out that of the particle description in the study of bulk 
properties of clusters and walks. The particle formalism becomes preponderant away from crit¬ 
icality, where conformal symmetry is lost but is replaced, in most of the interesting cases, by 
integrability, i.e. exact solvability of the particle dynamics. We saw that integrability is crucial 
in the non-perturbative study of near-critical singular behavior and of crossover phenomena. It 
plays a somewhat different role in problems, such as phase separation, in which the system is 
observed on scales much larger than the correlation length. In this case the order parameter 
profile and the fluctuations of the interface are determined exactly by general low-energy prop¬ 
erties of two-dimensional field theory, but in general the analysis requires some information on 
bound states, which is provided by integrability and determines the wetting properties, both in 
the bulk and at boundaries. 

We also exhibited the important role played by the fact that in two dimensions cluster bound¬ 
aries and interfaces can be associated to particle trajectories in imaginary time. Remarkably 
this leads, in particular, to a manifestation in real space of properties of the underlying particle 
description. So, the contact angle of phenomenological wetting theory corresponds to the posi¬ 
tion of a bound state pole, the dependence on the opening angle for wetting in a wedge (wedge 
covariance) reflects relativistic covariance, the asymptotic amplitude of crossing probabilities 
in near-critical percolation is determined by particle-boundary scattering. The fact that these 
and other results have been obtained in the last few years suggests that additional interesting 
developments can reasonably be expected for the future. 
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